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ABSTRACT 

Using very deep Spitzer MIPS 24 fim observations, we present an analysis of the bolometric lumi- 
nosities and UV extinction properties of more than 200 spectroscopically identified, optically selected 
(U n GlZ) z ~ 2 galaxies in the GOODS-N field. The large spectroscopic sample of rest-UV selected 
galaxies is supplemented with photometrically identified near-IR-selected ( "BzK" and "DRG" ) galax- 
ies and sub-mm sources at similar redshifts in the same field, providing a representative collection 
of relatively massive (M* > 10 10 M Q ) galaxies at high redshifts. We focus on the redshift range 
1.5 < z < 2.6, for which the 24 /im observations provide a direct measurement of the strength of the 
mid-IR PAH features in the galaxy spectra; the rest-frame 5 — 8.5 /im luminosities (Ls-s.s^m) are 
particularly tightly constrained for the objects in our sample with precise spectroscopic redshifts. We 
demonstrate, using stacked X-ray observations and a subset of galaxies with Ha measurements, that 
L5-s.5p.rn provides a reliable estimate of Lir, for most star forming galaxies at z ~ 2. We show that the 
range of Lir, in the optical/near IR-selected samples considered extends from ~ 10 10 L Q to > 10 12 L Q , 
with a mean (Lir) ~ 2 x 10 11 L©. The LIRG population at z ~ 2 is essentially the same population of 
galaxies that are selected by their optical/near-IR colors. Objects with LIRG to ULIRG luminosities 
are present over the full range of stellar masses in the samples, from 2 x 10 9 M Q to 5 x 10 11 Mq. We 
use the MIPS 24 /xm observations for an independent examination of dust extinction in high redshift 
galaxies, and demonstrate that, as in the local universe, the obscuration ( j~ IR ) is strongly dependent 
on bolometric luminosity, and ranges in value from < 1 to ~ 1000 within the sample considered. 
However, the obscuration is ~ 10 times smaller at a given Lboi (or, equivalently, a similar level of 
obscuration occurs at luminosities ~ 10 times larger) at z ~ 2 than at z ~ 0. We show that the 
values of Lir and obscuration inferred from the UV spectral slope (3 generally agree well with the 
values inferred from L5_g.5 Mm for Lboi < 10 12 L Q . As found previously by several investigators, for 
"ultraluminous" objects with Lboi > 10 12 L it is common for UV-based estimates to underpredict 
Lm by a factor of ~ 10 — 100. Using the specific star formation rate of galaxies (SFR per unit stellar 
mass) as a proxy for cold gas fraction, we find a wide range in the evolutionary state of galaxies at 
z ~ 2, from galaxies that have just begun to form stars to those which have already accumulated most 
of their stellar mass and are about to become, or already are, passively-evolving. 
Subject headings: cosmology: observations — dust, extinction — galaxies: evolution — galaxies: high 
redshift — galaxies: stellar content — infrared: galaxies 



1. INTRODUCTION 

The most direct method currently available for trac- 
ing the bolometric luminosities of high redshift star- 
forming galaxies (z > 2) has been from their submil- 
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limeter emission (e.g.. [ Smail et al1ll997t iHughes et alJ 
1998; Bar ger et al.lll998t) . Unfortunately, current sensi- 
tivity limits of bolometers and submillimeter wave in- 
terferometers allow for only the most luminous star- 
burst galaxies to be detected at high redshifts via their 
dust emission. Further compounding the problem is the 
coarse spatial resolution provided by such instruments, 
making it difficult to distinguish the counterpart (s) to 
the submillimeter emission for subsequent followup, al- 
though the recently developed method of radio-detection 
has been a breakthrough in alleviating this problem for 
most, but not all, b right submillimeter galaxies (e.g., 
Chap man et alJl2005() . Regardless, the dust properties 
of the vast majority of star-forming galaxies at high red- 
shift remained uninvestigated until recently. 

The rest-frame far-infrared (FIR) wavelength region 
is still inaccessible for the typical galaxy at redshifts 
z > 1, so we must look to other portions of the spectrum 
to directly examine dust properties. Our understand- 
ing of the mid-IR properties of local and high redshift 
galaxies advanced considerably with the launch of the 
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Infrared Space Observatory (ISO), which was sensitive 
enough to detect the mid-IR emission of 10 11 galax- 
ies at z ~ 1 (e.g-lFlores et al.lll999HElbaz et al J 120021: 
iPozzi et all 12004 iRowan-Robinson et al.ll2004|) . These 
observations revealed the almost ubiquitous presence of 
mid-IR dust emission features in star forming galaxies 
in both the local and z ~ 1 universe, and suggested the 
possibility of using the mid-IR dust emission of galaxies 
as a tracer of bolometric l uminosity (IBoselli et alJl!998t 
l Adelberger fe Steidelll2000l iDale et alJl2000HHelou et all 
boOOt lF5rster Schreiber et alJl2003|) 

These advances now continue with the highly success- 
ful Spitzer Space Telescope, providing the same sensi- 
tivity as ISO in probing dust emission from 10 11 L Q 
galaxies at z ~ 2. The progress made with Spitzer is 
particularly important for studying galaxies at z ~ 2 
because this epoch was until recently largely uninvesti- 
gated, yet is believed to be the most active in terms of 
star formati on and the build up of stellar and black hole 
mass (e.g.. IDick inson et al.M200 3t IRudn ick et al]|2003t 
Madau et al.|ll996HLillv et alJll996t ISteidel et1d]ll999l: 
Shaver et alJ 119961 iFan et all l200lt iDi Matteo et all 
2003HGiavalisco et al.lll996|) . The sensitivity afforded by 
the Spitzer MIPS instrument allows us to examine the 
typical L* galaxy at z ~ 2, rather than a limited slice 
of the most luminous population, a problem which, as 
alluded to before, limits the usefulness of submillimeter 
observations. 

It is fortuitous that the rest-frame mid-IR spectral 
features observed in local and z ~ 1 star-forming 
galaxies are redshifted into the Spitzer IRS spectral 
and MIPS imaging passbands at z ~ 2. The mid- 
IR spectral region from 3 — 15 (im is rich with emis- 
sion lines believed to arise from the stochastic heat- 
ing of small dust grains by UV photons (see review by 
iGenzel fc Cesarskvl l2000|) . These unidentified infrared 
bands (UIBs) are generally attributed to the C = C and 
C — H stretching and bending vibrational modes of a 
class of carbonaceous molecules called polycyclic aro- 
matic hydrocarbon s (PAHs; e.g., iPuget fc Legerlll989t 
iTielens et aT]ll999ft . which we assume hereafter. In the 
typical spectrum of a star-forming galaxy, these PAH 
emission lines, along with vario us fine-structure m etal 
and HI recombination lines (e.g., ISturm et alJ EOQOL are 
superposed on a mid-IR continuum thought to result 
from dust emission from very small grains, or VSGs 
ijDesert et aLl lr^fflO - ). In star- forming galaxies, the global 
PAH emission is mainly attributed to UV radiation 
from OB stars and has been f ound to correlate with 
global st ar formation rate (e.g., |Fjksl^-Sxhreib£r_el-aI. 
| 2004bl iForster Schreiber et alJ 12003 iRoussel et al 
2001), although variations with ioni zing intensity and 
meta lli city are also obse r ved (e.g.. lEngelbracht et al.l 
20051 IHogg et. all l20n3 IMonso-Herrero et al] 1200? 
Hclou et a,l.H200lUNorma,r,d et al.111 9951). 

Until now, the only way to estimate the bolometric 
luminosities of most galaxies at z ~ 2 independent of 
extinction was via their stacked X-ray and radio emis- 
sion: unfortunately these data are not sufficiently sen- 
sitive to detect individual L* galaxies at z ~ 2 (e.g., 
Rcdd v fc Steiddll200l INandra et 3,1112001 IBra,ndt et all 
12001ft . The Spitzer data considered in this paper are use- 
ful in assessing the bolometric luminosities of galaxies on 
an individual basis. One is still limited because detailed 



mid-IR spectroscopy is fea sible only for the most lumi- 
nous galaxies at z ~ 2 (e.g. JYan et al.l20T)5llHouck et alJ 
|2005|) . but L* galaxies at z ~ 2 (with L bo i ~ 10 11 L Q ) 
can be detected in deep 24 /xm images. We employ 
MIPS 24 /iin data to study the rest-frame mid-IR prop- 
erties of optical and near-IR selected galaxies at redshifts 
1-5 ^ z < 2.6. We describe the optical, near-IR, X-ray, 
and mid-IR data in § [21 and [3J Our large sample of spec- 
troscopic redshifts for optically-selected galaxies allows 
us to very accurately constrain the rest-frame mid-IR 
fluxes of z ~ 2 galaxies. In § we describe our method 
for estimating photometric redshifts for near-IR samples 
of galaxies where spectroscopy is less feasible. The pro- 
cedure for estimating infrared luminosities from MIPS 
data is outlined in § 15.11 We discuss the infrared lumi- 
nosity distributions of 24 /iin detected and undetected 
sources in § 15. 21 and 15. 31 The dust attenuation properties 
of optical and near-IR selected z ~ 2 galaxies and the 
correlation of these properties with bolometric luminos- 
ity are discussed in § and § The stellar populations 
and composite rest-frame UV spectral properties of faint 
24 /im galaxies are discussed in § [HJ In §[!|]we examine in 
more detail the mid-IR properties of massive galaxies at 
z ~ 2. We conclude in §EHby discussing the viability of 
optical and near-IR color criteria in selecting LIRGs and 
ULIRGs at z ~ 2 and what the Spitzer MIPS observa- 
tions can reveal about the mass assembly of galaxies at 
high redshift. A flat ACDM cosmology is assumed with 
H = 70 km s- 1 Mpc" 1 and fi A = 0.7. 

2. SAMPLE SELECTION AND ANCILLARY DATA 

2.1. Optical and Near-IR Selection 

The star-forming galaxies studied here were drawn 
from the sample of z ~ 2 galaxies in the GOODS- 
N field selected based on their observed U n G7Z col- 
ors to TZ = 25.5 l|Adelberger et all 12001 ISteidel et all 
12004]) . The optical images used for the selection of can- 
didates cover an area 11' by 15'. We refer to "BM" and 
"BX" galaxies as those which are selected to be at red- 
shifts 1.5 < z < 2 an d 2.0 < z < 2.6, respectively 
lAdelberger et al.ll200l ISteidel et all 12004]) . In addi- 
tion to BX/BM galaxies, we also c onsider galaxies se - 
lected using the z ~ 3 LBG criteria ijSteidel et alj f2003). 
The BX/BM and LBG candidates make up our U n G7Z, 
or optically, selected sample. We obtained rest-frame 
UV spectra with the blue channel of the Low Resolu- 
tion Imaging Spectrograph (LRIS-B) on Keck I for 386 
U n GTZ candidates. The numbers of candidates and spec- 
troscopically confirmed galaxies in the U n GTZ sample are 
summarized in Tabled The spectroscopic redshift distri- 
bution of U n GTZ galaxies in the GOODS-N field is shown 
in Figure F° r efficiency, we preferentially targeted for 
spectroscopy those U n GlZ candidates with 7?.-band mag- 
nitudes in the range TZ = 22.5 — 24.5 (AB units) and gave 
lower priorities for fainter objects where redshift identifi- 
cation is more difficult from absorption lines and brighter 
objects where the contamination fraction (from low red- 
shift interlopers) is larger. The star formation rate dis- 
tribution of spectroscopically-confirmed U n G7Z galaxies 
is similar to that of all U n G7Z g alaxies in the targ eted 
redsh ift range to K R = 21 (Vega: iReddv et alll2005t see 
also ISteidel et al.M200"4l for a discussion of spectroscopic 
bias of galaxy properties with respect to the photometric 
sample of U n G7Z galaxies). 
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TABLE 1 
Properties of the Samples 





















Sample 


Limits 


iV c a 




(arcmin 2 ) 










U n GK (BX/BM/LBG) 


72. < 25.5 


1571 


2.3 


10.2 ±0.3 


313 


2.25 ±0.33 


219 


0.65 


BzK /SF 


K s < 21.0 


221 


32 


3.1 ± 0.2 


53 


2.09 ±0.34 


82 


0.82 


Bzit/PE 


K B < 21.0 


17 


4 


0.24 ± 0.06 





1.70 ±0.20 


13 


0.54 


DRG h 


Ks < 21.0 


73 


19 


1.0 ±0.1 


5 


2.48 ±0.35 


24 


0.71 



a Number of candidates. 

b Number of directly-detected X-ray sources, including spectroscopically-confirmed galaxies. 

c Surface density of candidates. Errors are computed assuming Poisson statistics. 

d Number of spectroscopically confirmed objects with z > 1.5. Note that we only obtained spectra 
for those BzK and DRG galaxies that satisfy the U n G1Z criteria. 

c Mean redshift of sample for z > 1.5. For the BzK /PE sample, t his is the mean redshift of the 
photometric redshift distribution observed for BzK /PE galaxies (e.g. JDaddi et al.1 2004 Rcd dv et all 
2005). For the DRGs, this is the mean redshift of the spectroscopic redshift distribution observed 
for DRGs with z > 1.5 in four of the fields of the z ~ 2 optical survey [Rcd d7et~aT]|200g) . 

f Number of non-AGN galaxies (i.e., those with no direct X-ray detections) with spectroscopic 
redshifts 1.5 < z < 2.6. For the BzK /SF sample, this number includes both spectroscopically- 
confirmed BzK /SF galaxies (all of which are in the U n G1Z sample) and those with secure photomet- 
ric redshifts. For the BzK /PE sample, this includes all candidates without direct X-ray detections. 
For the DRG sample, this number includes galaxies with photometric redshifts 1.5 < z < 2.6. 

s Fraction of MIPS 24 (im detections to 8 /tjy (3 cr) among the 7V g galaxies. 

h Thc DRG samp le includes b oth star-forming galaxies and those with little star formation. As 
discussed in Rcddv et al. (2005), those with spectra (i.e., those DRGs which also satisfy the U n G1Z 
criteria) are likely to be currently forming stars. 



Our deep near-IR J and TV-band imaging, in ad- 
dition to publicly available deep optical imaging, al- 
lows us to select both star-forming galaxies and those 
with little current star formation in the GOODS-N 
field. Detai ls of the opt i cal an d near-IR images are 
provided in iReddv et all l|2005|) . The near-IR selec- 
ti on of star-forming galaxies is done using the criteria 
of lDaddi et all <l2f)04D. resulting in a "BzK /SF" sample 
(e.g.. lDaddi et al.l2004tlR.eddv et al.l2005jL The near-IR 
selected samples of galaxies with very little current star 
formation are constructed by considering the BzK and 
J — K s colors of candidates satisfying the BzK /PE and 
Distant Red Galaxy (DRG) criteria (e.g.. I Reddv et alJ 
120051 iDaddi et all 12004 iFranx et alJ 12003^ . Approxi- 
mately 70% of DRGs to K s = 21 (Vega) have signa- 
tures of intense star formation llPapovich et alJ 120051 
IReddv et al.l2005Uvan Dokkum et all2004|) *. The J and 
K band images cover a large area by near-IR standards 
(~ 8.5' x 8.5'), but are still less than half the area of the 
optical U n G7Z images. The number of near-IR selected 
candidates and their surface densities are presented in 
Tabled For the remaining analysis, we use AB units for 
optical (U n GlZ) magnitudes and Vega units for near-IR 
(K s ) magnitudes. 

2.2. X-Ray Data 

Th e very deep Chandra X -ray data in the GOODS-N 
field l) Alexander et alJ [2003) allow for an independent 
means of assessing the presence of AGN in the samples, 
which c an be quite significant for near-IR selected 
samples (Red dv et al.l l2005). In addition, we can stack 
the X-ray data for those galaxies lying below the X-ray 
detection thresh old to determine t h eir average emission 
properties (e.g. lLaird et al.l 200| ILehmer et all [2005: 

8 In order to ensure our sample is complete, we only consider 
DRGs to K B = 21. 




z 



Fig. 1. — Spectroscopic redshift distribution of optically (i.e., 
U n G1Z) selected galaxies, with typical error of cr(z) = 0.002. Also 
shown are arbitrarily normalized photometric redshift distributions 
of BzK /SF galaxies that do not satisfy the U n G7Z criteria, and 
DRGs. The errorbars indicate the average uncertainty in photo- 
metric redshifts (<x(z) = 0.45) for the near-IR selected BzK galax- 
ies and DRGs. 

I Reddv fc Steideill200l INandra et aDl200l iBrTndt et alJ 

2001). The X-ray dat a and stacking method s are 
discussed in deta il in IReddv fc Steidell {2004) and 
IReddv et alJ 1)20051) . The numbers of directly-detected 
X-ray sources in each sample considered here are sum- 
marized in Tabic n 

3. MID-IR DATA 

The mid-IR data are obtained from the Spitzer Multi- 
band Imaging Photometer (MIPS) instrument. The 24 
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micron data were taken as part of the GOODS Legacy 
Survey (P.I.: M. Dickinson) between May 27 and June 
6, 2004. They consist of 24 AORs (Astronomical Ob- 
servation Requests) of approximately 3 hours each. The 
combined data reach a depth equivalent to ~ 10 hours in- 
tegration at any point in the mosaicked image. The data 
are publicly available since August 2004 in the Spitzer 
archive. The basic calibrated data (BCD) produced by 
the Spitzer pipeline were used as the starting p oint for 
the d ata reduction. As explained in detail by iFaddal 
(2006), there were several artifacts which added noise 
to the images, hampering the detection of faint sources. 
These artifacts include image latencies from previous ob- 
servations of bright objects or image dark spots present 
on the pick-off mirror that are projected in different po- 
sitions by the cryogenic scan mirror during observations. 
Other variations come from the variable zodiacal light. 
We have corrected each BCD i mage for these effects us- 
ing the procedure described bv IFaddal ((2006) . The final 
mosaic consists of 7198 BC Ds combined using MOPEX 
(jMakovoz fc Marleaul l2005|) . The final reduced 24 /im 
mosaic of the GOODS-N region has a pixel size of 1'.'275 
and covers the entirety of our optical U n GlZ images and 
the measured 3 a depth is approximately 8 /(Jy. The 
large beam-size of MIPS (~ 5'.'4 at 24 /im) combined 
with the effects of blending make aperture photometry 
impractical. Instead, we have chosen to use a PSF fitting 
method to extract 24 /im fluxes for our galaxies, similar 
to the method used to extract IRA C fluxes for galaxies 
as discussed in lShaplev et al.1 l|2005l h 

We first compute an empirical PSF using a two-pass 
approach. In the first pass, we take the median flux pro- 
file of several tens of distinct MIPS point sources across 
the mosaic to create a rough empirical PSF. This PSF is 
then used to subtract sources in proximity to the tens of 
MIPS point sources in order to better determine the PSF 
profile out to larger radii. This is important for the MIPS 
PSF which contains up to 20% of the point source flux in 
the Airy rings. The effects of source confusion are miti- 
gated by employing the higher spatial resolution Spitzer 
IRAC data in the GOODS-N field to constrain the MIPS 
source positions. The empirical PSF, normalized to unit 
flux, is fit to these positions and the fluxes are extracted. 
The PSF used here extends to 15'.'3 radius, and we apply 
a 15% aperture correction based on the observed curves- 
of-growth of MIPS point source profiles from the First 
Look Survey. 

The number of MIPS detections (> 3 cr) and non- 
detections in each sample are summarized in Table 
Virtually all of the directly-detected X-ray sources are 
detected at 24 /mi and all have optical/X-r ay flux ratios 
indic ating that AGN dominate the X-flux l)Reddv et alJ 
|2QQ5t). Submillimeter galaxies are often associated 
with direct X-ray d etections (e.g.. I Alexander et alJl2005l 
IReddv et all [2005) even though their bolometric lumi- 
nosities may still be dominated by star-formation. Since 
we are primarily interested in the mid-IR emission as a 
tracer of star formation, we have excluded all directly- 
detected X-ray so urces (almost all of which are AGN; 
IReddv et al.ll2005j) for most of the analysis, unless they 
happen to coinc i de wit h a radio-detected SMG from the 
iChaoman et all 1)20051) (SC05) sample. We caution that 
the resulting sample of 9 SMGs is not meant to be 
uniform or complete: ~ 40% of SMGs are not associ- 



TABLE 2 
Local Template Galaxies 



Name a 


z h 


(xlO 10 L ) 


M82 


0.000677 


1.8 


NGC253 


0.000811 


1.8 


M83 


0.001711 


1.9 


M51 


0.002000 


2.1 


NGC1808 


0.003319 


3.8 


NGC1097 


0.004240 


3.8 


NGC1365 


0.005457 


8.7 


NGC520 


0.007609 


6.5 


NGC7714 


0.009333 


5.6 


NGC3256 


0.009386 


40 


NGC6240 


0.024480 


60 


IRS9 


1.83 


1800 



a All galaxies are compiled from 
I Forster Schreiber et al.l 12004al) and 
iRoussel et alJ 120011) . except for IRS9 
and NGC7714 w hich are taken from 
lYan et alJ 120051) . The bolometric 
luminosity of IRS9 is c onstrained 
within a factor of 2-3 lYa n et al.1 

b Redshifts for local galaxies arc 
obtained from NASA/IPAC Extra- 
galactic Database (NED). 

ated with radio sources, either b ecause of their higher 
redshifts or radio faintness (e.g., IChaoman et alJ 120051: 
ISmail et al.M2002|) . Furthermore, the submillimeter ob- 
servations are not uniform over the field. Nonetheless, 
it is of obvious interest to at least compare the mid-IR 
properties of this limited set of radio-detected SMGs to 
those of galaxies in other samples. 

The MIPS 24 /im filter directly samples the rest-frame 
luminosity from 5 — 8.5 /im (Ls-g.s^m) for redshifts 
1-5 < z < 2.6. We used the mid-IR spectral shapes of 
star-forming galaxies (listed in Table |5J) as templates 
in order to if-correct the 24 /tm fluxes to determine 
£5-8.5^111- Figure El shows the expected 24 /tm fluxes of 
the galaxies listed in Table [21 as a function of redshift. 
These galaxies were chosen to cover a large range in 
SFRs (from quiescent spiral galaxies to starbursts and 
LIRGs/ULIRGs). The mid-IR s pectra are obtained from 
either ISO (for local galaxies; IForster Schreiber et alJ 
l2004al) or Svitzer IR S (for z ~ 2 hyper-luminous galax- 
ies: lYan et al.ll2005fl observations. Some properties of 
these galaxies are summarized in Table The mid-IR 
spectrum of each galaxy is redshifted, convolved with 
the MIPS 24 /tm filter, and the fluxes are normalized to 
have i 5 _ 8 .5 Mm = 10 10 L to produce the dotted curves 
111 Figure □ The small dispersion between the templates 
over redshifts 1.5 < z < 2.6 reflects small changes in 
the broadband 24 /im fluxes of galaxies due to changes 
in the relative strengths of the various PAH emission 
lines and the ratio of PAH-to-continuum flux. The solid 
curve in the figure shows the average of all the template 
galaxies. The prominent peak at z s=s 1.9 is primarily due 
to the 7.7 /tm feature shifting into the MIPS 24 /tm filter. 

4. PHOTOMETRIC REDSHIFTS OF NEAR-IR SELECTED 
GALAXIES 
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FlG. 2. — Expected 24 fim fluxes of the local and high redshift template galaxies summarized in Table 121 shown by the dotted curves, 
normalized so that Ls—s.Sfim — 10 10 Lq. The solid curve is the average over all the templates. Also shown are the observed fluxes of 
U n G1Z (small points) and BzK/SF galaxies (thin triangles) with spectroscopic redshifts; and BzK/SF (thick triangles), BzK/PF, (solid 
pentagons), and DRGs (open squares) with photometric redshifts. The typical err or in photometric reds hifts of DRG and BzK galaxies 
is a(z) ~ 0.45. Radio-detected submillimeter galaxies to Sgso^m ~ 4.2 mjy from Chapman ct al. (2005) are shown by the open circles. 
We have removed directly-detected X-ray sources in the samples above, except for those in the SMG sample. Crosses denote hard-band 
X-ray sources with spectroscopic redshifts in the U n GlZ sample. The vertical errorbar in the upper left-hand side of the figure shows the 
typical uncertainty in the 24 /an flux. The shaded region indicates the redshift range over which the MIPS 24 (im filter directly samples 
the rest-frame 7.7 /an PAH feature. 



Figure [3 illustrates the sensitive dependence of the K- 
correction on the redshift (e.g., galaxies with a given ob- 
served /24^ m can have a factor of 3 spread in £5-8.5^111 
depending on their redshift in the range 1.5 < z < 2.6). 
Our large spectroscopic sample gives us the distinct ad- 
vantage of knowing the precise redshifts for the optically 
selected galaxies, removing the added uncertainty intro- 
duced by photometric redshifts where the precise loca- 
tion of the PAH features with respect to the MIPS 24 fim 
filter is unknown, adding consider able uncertainty to th e 
inferred infrared luminosities fe.g,. IPapovich et all 2005). 
As we show below, the typical error in the photometric 
redshifts derived for near-IR selected galaxies (even when 
using data across a large baseline in wavelength, from UV 
through Spitzer IRAC), is o~(z) ~ 0.5. This error in red- 
shift translates to at least a factor of three uncertainty 



in Lk_ 



5— 8.5/x m- 



Nonetheless, photometric redshifts are the only practi- 
cal option for optically faint galaxies where spectroscopy 
is not feasible. This is true for many of the DRGs 
and BzK/VEj galaxies. We supplement our spectro- 
scopic database of optically selected galaxies with pho- 
tometric redshifts of near-IR selected galaxies. We made 
use of the HyperZ code to d etermine photometric red- 
shifts (Bolzon ella et alJl2000|) . To quantify the uncer- 
tainties in photometric redshifts, we tested the code on 
U n GlZ galaxies with spectroscopic redshifts, fitting to the 
U n BGVTlIz + JK photometry. The BVIz magnitudes 
are obtai ned from the vl.l relea se of the GOODS ACS 
catalogs ijGiavalisco et al.| [2"004). Errors in the optical 
U n GTZ and near-IR JK magnitudes a re determined from 
Mo nte Carlo simulations described bv lErb et alJ ((2006b) 
and lShaplev et alJ l)2005|) . 

The x 2 between the modeled and observed colors was 
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Fig. 3. — Left panel: Photometric versus spectroscopic redshifts for U n GTZ galaxies. The solid line indicates z p hot = z Bpcc . Right panel: 
Az, as defined in the text, as a function of stellar mass for galaxies with redshifts 1.5 < z < 2.0 (solid circles) and 2.0 < z < 2.6 (open 
circles). The figure only shows objects with > 3 a detections in at least two IRAC channels. 



calculated for each galaxy for a number of different star 
formation histories (with exponential decay times r = 10. 
20, 50, 100, 200, 500, 1000, 2000, and 5000 Myr, and 
t = oo) as a function of redshift. The best-fit photomet- 
ric redshift is the redshift at which \ 2 is minimized. As 
a figure of merit of the resulting fit (and for easy com- 
parison with other studies), we compute 

Az = (Zphot — Z spcc )/(1 + Zgpoc). (1) 

The results are shown in the left panel of Figure |3] The 
dispersion in Az is a(Az) « 0.15 for galaxies with spec- 
troscopic redshifts 1.5 < z < 2.6. The actual error in 
redshift is typically a(z) ~ 0.45. Both thelShapley et all 
l)2005[) code ( which us e slBruzual fc Chariot! (|2003f) mod- 
els) and the IBenfted l|2000fl code gave similar results. 
We found that adding Spitzer IRAC data does little to 
tighten the scatter between photometric and spectro- 
scopic redshifts for most galaxies, reflecting the absence 
of any distinct features (e.g., strong spectral breaks) 
across the IRAC bands for z ~ 2 galaxies. We note 
from the left panel of Figure that photometric red- 
shifts systematically underestimate the true redshifts of 
z > 3 galaxies. This should not significantly affect our 
subsequent analysis since we only consider galaxies up 
to z = 2.6, and most of the BzK/SF and DRG galaxies 
have photometric redshifts z p hot ^ 2.5. 

The IRAC data are nonetheless a powerful tool in dis- 
cerning the more massive galaxies from the less mas- 
sive ones. Since most of the optically- faint DRGs and 
BzK /PE galaxies are on average amongst t he more mas- 
sive g alaxies at z ~ 2 (e.g., Figure 18 of iReddv et al J 
( 2005)), we have incorporated the IRAC data in the pho- 
tometric redshift fits. The right panel of Figure [21 shows 
Az as a function of inferred stellar mass for U n GTZ galax- 
ies computed for the best-fit r model, normalizing to the 
observed colors 9 . The scatter in Az for galaxies with stel- 
lar masses M* > 10 10 M Q is er(Az) ~ 0.20, and we will 
assume this value for the error in photometric redshifts of 

9 We assume a Salpeter IMF in calculating the stellar mass. 



the DRG and B zK/PE galaxies. To extend the compari- 
son presented bv IReddv et al.1 l|2005|) between U n GlZ and 
BzK selected samples of star-forming galaxies by exam- 
ining their mid-IR properties, we compute photometric 
redshifts for BzK /SF galaxies that do not satisfy the 
U n GlZ criteria. For the BzK /SF galaxies, we assume an 
error of <j(Az) ~ 0.15, according to Figure^. 

We obtained 51 secure photometric redshift fits for 
BzK /SF galaxies not selected by the U n GlZ criteria (out 
of 95 such objects), and their (arbitrarily normalized) 
photometric redshift distribution is shown in Figure ^ 
Also shown is the (arbitrarily normalized) photometric 
redshift distribution for 28 (out of 49) non X-ray de- 
tected DRGs for which we were able to derive secure 
photometric redshifts 1 < z p i lot < 3.5. The remaining 
DRGs either have z p hot < 1 or z p hot > 3.5 (and are irrel- 
evant to the analysis considered bel ow) or had photom- 
etry t hat was inconsistent with the iBruzual fc Chariot! 
( 2003) models considered here, resulting in a large x 2 
value between the model and observed colors. The DRGs 
examined here appear to span a very large range in red- 
shift from z ~ 1 — 3.5 , a result similar to that found by 
IPanovicn et, piol for DRGs in the GOODS-South 
field. We obtained good photometric redshift fits for 
only two of the BzK/PE galaxies: one at z ~ 1.2 and 
the other at z ~ 1.5. We reiterate that for purposes 
of redshift identification, we only assumed the photo- 
metric redshifts for those galaxies for which we were 
able to obtain good SED fits (i.e., with x 2 w 1) to 
the observed data. There were a number of objects 
for which the photometric redshift errors were relatively 
large (Sz/(1 +z) > 1) or had derived redshifts that were 
much larger (z > 4) or smaller (z < 1) than of interest 
here, and we excluded such objects from our analysis. 
Hence, the photometric redshift distributions in Figure^ 
should not be attributed to the populations as a whole. 
For the remaining analysis we consider only optically- 
selected galaxies with spectroscopic redshifts and near- 
IR selected galaxies with photometric redshifts between 
1-5 ^ z < 2.6 where the 24 fim fluxes directly trace the 
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flux at rest-frame 7.7 /im. This is indicated by the shaded 
region in Figure 

5. INFRARED LUMINOSITIES OF OPTICAL, NEAR-IR, 
AND SUBMILLIMETER SELECTED GALAXIES AT Z ~ 2 

5.1. Inferring Infrared Luminosities from Ls-g.S/jm 

The conversion between Ls-g.Suni and infrared lumi- 
nosity will largely depend on the assumed spectral tem- 
plate relating the mid-IR emission of galaxies to their 
total infrared luminosities. Fortunately, the deep X-ray 
data allow us to determine whether L5_s.5/im scales with 
infrared luminosity (or star formation rate) independent 
of any assumed template, adopting the local empirical 
relationship between X-ray and FIR luminosity for star- 
forming galaxies. Figure 01 shows the ratio of I<5-8.5^m 
to stacked X-ray luminosity of (X-ray undetected) galax- 
ies in bins of £5-8.5^111 : we only considered optically- 
selected galaxies with spectroscopic redshifts since it is 
for these galaxies which we are able to most accurately 
constrain the rest-frame X-ray luminosities. Since X-ray 
emission is sensitive to star formation on time-scales of 
> 100 Myr (see §EJ, Figure 01 shows results excluding 
galaxies with inferred ages < 100 Myr. Each bin contains 
10 — 20 sources except the faintest bin which includes 
45 galaxies undetected at 24 /im with ages > 100 Myr. 
The X-ray data for galaxies in each bin were stack ed us- 
ing the procedure described in lReddv et al.1 ((2005). The 
mean value of the mid-IR-to- X-ray luminosity ratio is 
(L5-S. Bfua/L^-io kcv) ~ 251 ± 41. The X-ray luminosi- 
ties of local star-forming galaxies are found to tightly 
correlate with their infrared emission for galax ies with 

10 8 < L Fm < 10 12 L fe.g.. lRatia11i7t^l200l . Using 
the X-ray luminosity as a proxy for infrared luminos- 
ity therefore implies that the rest-frame mid-IR fluxes 
follow the total infrared luminosity (Lir) for the vast 
majority of optically-selected galaxies at z ~ 2 10 . The 
mid-IR fluxes must also follow the infrared luminosity 
for most near-IR selected star-forming galaxies as well 
given the large overlap (70% — 80%) between optical and 
near-IR selected sa mples of z ~ 2 star-forming galaxies 
(IReddv et al.l l2005'l. As we show below, the conversion 
we assume between rest-frame mid-IR and infrared lu- 
minosities reproduces the average infrared luminosities 
predicted from stacked X-ray analyses (§ IB.lfl , 

To quantify the relationship between rest-frame mid- 
IR and total infrared luminosity with a minimum num- 
ber of assu mptions, we have made use of the data 
compiled by lElbaz et~ O <|2002), which include IRAS 
60 and 100 /an measurements and ISO observations of 
149 local star-forming galaxies with Lir in the range 

10 9 < iiR < 10 126 Li©. The mean and dispersion of 
the IR/MIR flux ratio for the sample of 149 galaxies 
is (log(Li R /i 5 -8.5 M m)) - 1-24 ± 0.35. Note the large 
1 a dispersion of a factor of 2.2 in the IR/MIR flux ra- 
tio. The dispersion in the IR/MIR flux ratios between 
galaxies may be driven partly by changes in the mid- 
IR line strengths as the aromatic carriers are dehydro- 
genated and/or destroyed dependin g on the intensity of 
the ambient UV ionizing field (e.g., EI 

onso-Herrero et all 

10 Another commonly used definition of Lir is the total luminos- 
ity from 1 — 1000 fim. This will differ from Lg_ 1000^™ by only a few 
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Fig. 4. — Ratio of L5— g.s^m and stacked rest-frame 2 — 10 keV X- 
ray luminosity in bins of L5— s.S/im for galaxies with redshifts 1.5 < 
z < 2.6 and those with inferred ages greater than 100 Myr. We 
have excluded sources directly-detected in the Chandra 2 Ms data. 
The faintest bin is for galaxies undetected at 24 fim. Horizontal 
error bars indicate the 1 <r dispersion in Ls_ a.Bum m each bin. 
Vertical error bars show the uncertainty in the mid-IR-to- X-ray 
luminosity ratio, computed as the uncertainty in the mean X-ray 
luminosity added in quadrature with the Poisson error in the mean 
mid-IR luminosity of galaxies in each bin. The dashed horizontal 
line and shaded region indicate the mean and 1 cr uncertainty of 
the mid-IR-to-X-ray luminosity ratio of ~ 251 ± 41. 
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percent, and for the remaining analysis, we take Lir 
as defined by Sanders & Mirabel ( 1996' 
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1995). Metallicity effects and a changing distribution of 
dust with respect to HII regions also likely contribute to 
the observed factor of 2 — 3 dispersion in the IR/MIR ra- 
tios. Nonetheless, the mean IR/MIR flux ratio is similar 
to that inferred from the lDale et al.l l)2001|) template SED 
for a median log[/iy(60^m)//, y (100^ m)1 ~ —0.20 . Base d 
on the sample of 149 galaxies from lElbaz et alJ {2002), 
we assume L m w 17.2L 5 _ 8 . 5Mm to convert L 5 - 8 . 5flm to 

It is worth noting that the relationship between 
-^5-8.5^111 and total IR luminosity for local star-forming 
galaxies may be described by a more complicated func - 
tion, such as a two power-law fit (e.g.. lElbaz et al.l 20021. 
These complicated relationships may not apply to galax- 
ies at z ~ 2 for several reasons. First, the IR/MIR ratio 
may change depending on the relative contribution of 
older stellar populations in heating PAH and larger dust 
grains. The heating of dust by cooler stars is expected 
to be more prevalent for the less luminous local galaxies 
with older stellar populations, on average, than for z ~ 2 
galaxies with relatively younger stellar populations. Sec- 
ond, it is found that z ~ 2 galaxies have metallicities that 
are on average 0.3 dex lo wer than those of local galax- 
ies at a given stellar mass (|Erb et al.ll2006a() . Therefore, 
the metallicity depen dence of the IR/MIR ra tio found 
for local galaxies (e.g-. lEngelbracht et alJl2005^ suggests 
that we may not be able to ascribe the IR/MIR ratio 
for a galaxy of a given stellar mass at z = to a galaxy 
at z ~ 2 with the same stellar mass. A third possibil- 
ity, and one that is suggest ed by the results of this pa- 
per and other studies (e.g ., lAdelberger fc Steidell I200fl 
iCalzetti &: Heckmanlll999|) . is that the dust obscuration 
of galaxies at a given bolometric luminosity changes as 
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function of redshift, a result that may reflect dust enrich- 
ment and/or a changing configuration of dust as galax- 
ies age. Therefore, the relative distribution of PAH and 
larger dust grains within galaxies may also change as a 
function of redshift. Because of these uncertainties, and 
since the primary motivation of our study is to indepen- 
dently establish the validity of MIPS observations to infer 
the infrared luminosities of z ~ 2 galaxies, we adopted 
the simplest conversion that assumes only that the typ- 
ical IR/MIR luminosity ratio for local galaxies with a 
wide range in infrared luminosity applies at z ~ 2. By 
taking an average over the local galaxies, we ensure that 
the derived Lir is not more than a factor of 2 — 3 away 
from that predicted using the IR/MIR ratio of any in- 
dividual galaxy. As we show below, our constant con- 
version reproduces within the uncertainties the results 
that we obtain from stacked X-ray and dust-corrected 
UV estimates of Lir. 

In addition to the stacked X-ray and dust-corrected 
UV estimates, we also have spectroscopic Ha observa- 
tions for a small sa mple of 10 U U GT Z galaxies in the 
GOODS-North field (|Erb et al. 2006c) with clean (i.e., 
not blended) MIPS detections. Once corrected for ex- 
tinction, the Ha fluxes of these galaxies provide estimates 
of their total (bolometric) luminosities, which we take to 
be the sum of the Lir and observed 1600 A luminosity 
(uncorrected for extinction): 

Lho\ = Lir + ii6oo- (2) 

In Figure [5] we show the comparison between Lboi esti- 
mated from the sum of the MlPS-inferred Lir and ob- 
served 1600 A luminosity (^boi'" UV ) w ^ n ^ b °i estimated 
from the spectroscopic Ha observations (Lj^j). The re- 
sults indicate that within the uncertainties the two esti- 
mates of iboi track each other very well (with a scatter 
of 0.2 dex) over the range of Lboi typical of galaxies in 
optical/near-IR selected samples (§01, at least for this 
limited sample of 10 objects. The agreement between the 
MIPS and Ha-inferred estimates suggests that our con- 
version relation between L5-8.5^tm and Lir works reason- 
ably well. Nonetheless, for comparison with our constant 
conversion relation, we also consider the effect on our re- 
sul ts of assuming a tw o power-law conversion suggested 
bv lElbaz et all {2002 1 As we show below, assuming the 
two power-law conversion does not change the main con- 
clusions of our study. 

The far-infrared luminosity (Lfir) is typi cally defined 
to be the luminosi t y fro m 40 — 120 /xm ijHelou et al.l 
dH). ISoifer et all Ip8l found L m ~ (1.91 ± 0.17) x 
Lfir for galaxies in their Bright Galaxy Sample. Model- 
ing of the warm and cool components of the dust emis- 
sion in U V-bright galaxies ind icates a conversion factor 
of ~ 1.75 l)Calzetti et al.ll2000j) . We take a median value 
of ~ 1.80 in converting the inferred Lir of galaxies to a 
FIR luminosity. Generally, uncertainties in the conver- 
sion between Lir and Lfir are much smaller than the 
uncertainties in converting Ls-g.s^m to Lir. 

Hereafter we assume uncertainties in the total infrared 
luminosities as follows. For U n GlZ galaxies and radio- 
selected SMGs with spectroscopic redshifts, we assume 
an uncertainty in log Lir of 0.3 dex, corresponding to 
the dispersion in the conversion between m and 

Lir. For near-IR selected BzK galaxies and DRGs, the 
photometric redshift error will add an additional 0.5 dex 
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Fig. 5. — Comparison between MlPS-inferred £™[ +UV and Ha- 
inferred for a sample of 10 UnGlZ galaxies. The dispersion 
in the ratio of the two estimates is ~ 0.2 dex for the subsample 
considered here. The errorbar shows the typical uncertainty of each 
point. 

scatter. The total uncertainty in log Lir for BzK galax- 
ies and DRGs with photometric redshifts is 0.6 dex. 

5.2. Infrared Luminosity Distributions 

Figure [2] summarizes the observed /24pm fluxes of 
galaxies as a function of redshift. In this figure, all 
direct X-ray detections were removed from the U n GlZ, 
BzK and DRG samples, unless they happened to coin- 
cide with a radio-det e cted s ubmillimeter galaxy (SMG) 
from iChapman et alJ ((2005) , or unless they have spec- 
troscopic redshifts in the U n GlZ sample (crosses in Fig- 
ure 01. The U n GlZ and SMG samples include objects 
outside the region covered by our near-IR imaging. The 
overlap betwe en the samples cons idered here is discussed 
extensively in iReddv et al.l l)2005[) . The observed 24 /an 
fluxes for objects in the various samples generally span 
a large range, from our sensitivity limit of ~ 8 /iJy to 
> 200 /iJy. 

For a more meaningful comparison, we have computed 
Lir for galaxies using the prescription described in § 15.11 
Figure shows the distribution of Lir as inferred from 
L5_8.5/im for galaxies in t he U„G1Z, BzK/SF, DRG, 
and radio-detected SMG ijChapman et alJ [2005) sam- 
ples with either spectroscopic or photometric redshifts 
1.5 < z < 2.6. As in Figure [3 we have excluded directly 
detected X-ray sources from the distributions unless they 
coincide with a radio-detected SMG source. We show the 
resulting distributions assuming the constant conversion 
and two power-law conversion relations in the left and 
right panels, respectively. The distributions assuming 
the two power-law conversion are bimodal. However, the 
distributions of observed UV luminosities and dust cor- 
rection fac tors of U„G1Z galaxies are approximately gaus- 
sian (e.g., ISteidel et xUMM IReHdv et nlJ 12(H)(1). As- 
suming the Calzctti ct al. ( 2000) law to convert the ob- 
served luminosities to dust-corrected luminosities then 
implies that the bolometric luminosity distribution of 
U n GlZ galaxies should be gaussian, a result not in accor- 
dance with the bimodal distribution computed assuming 
the two power-law conversion. More generally, we expect 
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Fig. 6. — Distribu tions of Lu as inferred from L5— g.s^m for galaxies in the U n G1Z, BzK/SF, DRG, and radio-detected SMG 
IChapm an ct al. 2005) samples with rcds hifts 1.5 < z < 2.6 assuming our constant conversion between L5-g.5 Mm and Lir (left panel) 
and the two power-law conversion of[Elbaz ct al. (2002) (right panel). We have excluded directly-detected X-ray sources unless they hap- 
pen to coincide with a radio-detected SMG. 



to find gaussian distributions of luminosity for galaxies 
in photometric surveys since such galaxies are typically 
selected on a continuous range of observed colors and / or 
magnitudes. The bimodality in Figure EJd is likely due 
to the sparse data used to establish the two power-law 
relation for galaxi es with 1 x 10 11 < Lir sa< 5 x 10 11 L Q 
(see Figure 5d of lElbaz et al.ll2002() . This range of L m 
happens to encompass the typical IR luminosity of z ~ 2 
galaxies as inferred from X-ray a nd dust-corrected UV es- 
timates l|Reddv fc Steide]|l2004j) and it is partly for this 
reason that we favored our simple constant conversion 
relationship. 

Regardless of the conversion used, we find that the 
bulk of the galaxies in the U n GH and BzK /SF samples 
and detected at 24 /jm have inferred infrared luminosities 
comparable to those of local luminous infrared galaxies 
(LIRGs), with 10 11 < L m < 10 12 L . Galaxies in the 
U n GlZ sample with fa^m si 8 /iJy (corresponding to the 
3 a sensitivity limit) have (Lir) ~ 3.1 x 10 11 L Q for the 
constant conversion and (Lir) ~ 2.1 x 10 11 L Q for the two 
power-law conversion (the two power-law distribution is 
broader than that obtained using the constant conver- 
sion). Both the U n GlZ and BzK /SF samples also host 
galaxies which, based on their inferred Lir, would be con- 
sidered ultra-luminous infrared galaxies (ULIRGs) with 
-^IR ^ 10 12 Lq. Note that if we excluded all direct X- 
ray detections, including the submillimeter sources, the 
maximum inferred Lir of U n GlZ and BzK /SF galaxies 
is w 2.5 x 10 12 Lq, an infrared luminosity which is simi- 
lar to the detection limit of th e Chandra 2 Ms data for a 
galaxy at z ~ 2 assuming the iRanalh" eta!] l)2003T) con- 
version between X-ray and FIR luminosity. 

The BzK /SF sample distribution shown in FigureElin- 
cludes galaxies that do not satisfy the U n GlZ criteria (i.e., 
BzK/SF-U n G1l galaxies). These galaxies (to K s = 21) 
have a mean IR luminosity that is identical to that of 
U n G7Z galaxies to K s — 21. The average IR luminosity of 
U n G7Z galaxies is ~ 1.8 times fainter than BzK galaxies 



since the U n GlZ sample includes galaxies which extend 
to fainter K s magnitudes than those in the BzK sample. 
Therefore, while the BzK /SF-U n G1Z galaxies have red- 
der G — 1Z colors than required to satisfy the U n GlZ cri- 
teria, it appears that their infrared luminosities are still 
comparable to those of U n G1Z galaxies (see also the dis- 
cussion in §EJ, a result consistent with tha t obtained in 
X-ray stacking analyses ijReddv et al.l2005l) . FigureElin- 
dicates the DRG galaxies with photometric redshifts be- 
tween 1.5 < z < 2.6 also span a large range in Lir, from 
luminosities characteristic of LIRGs to ULIRGs. We find 
a luminosity distribution of DRGs to K s = 21 that is 
in good agreement w ith the Lir distribution found by 
iPapovich et alJ l)2005|) for DRGs (to approximately the 
same depth) in the GOODS-South field 11 . We note that 
U n GlZ galaxies and DRGs to K s = 20 have the same 
inferred Lir as K s < 20 galaxies selected in other ways 
(e.g., using the BzK /SF criteria). 

The inferred Lir for the one BzK /PE selected galaxy 
with z ~ 1.5 is ~ 1.2 x 10 12 L Q . The mean /24 Mm flux 
of MlPS-detected (and non X-ray detected) BzK/VF 
galaxies without redshifts is {fiApm) ~ 81.4 / xJy which, at 
the mean redshift of BzK /VF galaxies fe.g. JReddv et all 
120051 iDaddi et ail l2004|) of z ~ 1.7, corresponds to 
Lir ~ 6 x 10 11 L . The 24 /an detection rate (~ 50%; 
Tabled of non X-ray detected BzK/YFi galaxies im- 
plies some contamination by star-forming galaxies; this 
is not unexpected given that photometric scattering can 
have a significant eff ect on samples con structed by color 
selection techniques l|Reddv et al.ll2005]) . 

The radio-detected submillimeter galaxies to SgsoAim ~ 
4.2 mJy analyzed here have inferred Lir of 10 12 < Lir < 
10 13 L Q , which can be directly compared with their 
bolometric luminosities ca lculated based on the s ubmil- 
limeter data presented by IChapman et alJ (J2005) . The 

11 The DRG sample of lPapovich et alJ <2005l) extends to K s = 
23.2 in AB magnitudes, or K B ~ 21.4 in Vega magnitudes, over an 
area twice as large as studied here. 
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850 /xm-inferred bolometric luminosities (L IR Mm ) are 
sensitive to the assumed characteristic dust temperature 
associated with a greybody fit to the dust SED and the 
emissivity. For example, a change in the assumed dust 
temperature from = 36 K (th e median value for the 
sample of radio-detected SMGs in lChapman et aI1l2005[) 
to a cooler temperature of Td = 31 K results in a factor 
of ~ 5 decrease in the inferred FIR luminosities. Fig- 
ured shows the comparison between 850 /zm and 24 /im 
inferred bolometric luminosities of radio-detected SMGs. 
We also show th e point corresponding to IRS9 from the 
lYan et af] l)2005f ) sample of hyperluminous z = 2 sources 
with IRS spectroscopy — this source has independent 
constraints on its FIR luminosity based on MIPS 70 and 
160 /xm and MAMBO millimeter measurements. The 
infrared luminosity of IRS9 is Lir ~ 1.8 x 10 13 L Q (con- 
strained to within a factor of 2 — 3) based on these multi- 
wavelength measurements l)Yan et al.ll2005|) . 
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Fig. 7. — Lj^° Mm , inferred from the submillimeter fluxes of radio- 
detected SMGs IChapman et al1l2005l) . plotted against j^ 8 - 5 /" 71 
inferred from Ls_ s.5 pm assuming a IR/MIR flux c onversion of 
17.2 (s olid circles) and the two power-law conversion of lElbaz et alj 
(2002) (crosses). The infrared luminosity of IRS9 is inferred from 
MIPS 24, 70, and l&O^m data and MAMBO millimeter measure- 
ments lYan et alj 1200a) . Also shown are the local star- forming 
templates listed in Table |21 

Figure □ shows that the IR/MIR flux ratio for IRS9 is 
comparable (to within ~ 1 a) to those of the local star- 
forming galaxies listed in Tabled — these local galaxies 
are 1 — 3 orders of magnitude fainter than IRS9. Judg- 
ing the validity of our conversion for the hyper-luminous 
galaxies at z ~ 2 is difficult given that very few of these 
galaxies have independent constraints on their bolomet- 
ric luminosities. On the other hand, the submillimeter 
(850 /im) inferred infrared luminosities of bright SMGs 
are systematically a factor of 2 — 10 higher than pre- 
dicted using our conversion between Ls-g.s^m and Lir. 
The crosses in Figure [7| demonstrate that the systematic 
offset cannot be completely accounted for if we assume 
a two power-law conversion between Ls-s.s^m and Lir 
- there are still 4 of 9 SMGs with L^ m that are fac- 
tor of 2 — 10 larger than predicted from their 24 /jm 
fluxes, and the distribution of SMG points when con- 



sidering the two-power law conversion is not symmetric 
about the line of equality (solid line in Figure Q. One 
possibility is that the submillimeter estimates are in fact 
correct and that our assumed conversion between mid- 
IR and IR luminosities does not apply to SMGs. The 
IRS sample considered here consists of just one hyper- 
luminous galaxy at z = 2, and if we ignore this galaxy 
then the systematic offset of SMGs may indicate a break- 
down of our assumed conversion for the most luminous 
sources at redshifts z ~ 2 with Lir > 10 13 L Q . The sec- 
ond possibility is that the submillimeter estimates over- 
predict the infrared luminosities of SMGs and that our 
MIR-to-IR conversion is correct. This may not be sur- 
prising since the conversion between submillimeter flux 
and bolometric luminosity is very sensitive to the as- 
sumed dust temperature, and a decrease in the assumed 
temperature of just a few degrees can reduce the in- 
ferred bolometric luminosity by a factor of ~ 5 — 10 
(see example above). Finally, it is possible that nei- 
ther the submillimeter or mid-IR inferred infrared lu- 
minosities of bright SMGs is correct. We note that it is 
common for these luminous galaxies to host AGN, and 
this can alter the observed mid-IR and IR fluxes beyond 
what would be exp ec ted given pure star f ormation (e.g., 
[ Armus et ail 120041 iFadda et all l2002t lAlmaini et all 
1999; [Fabian k. Iwasawal I1999J) As another example, 
Arp 220 has anomalously low PAH emission for its bolo- 
metric luminosity (when compared with other ULIRGs), 
suggesting that the galaxy contains a buried quasar 
and/or a heavily dust-enshrouded starburst such that 
the extinction at rest-frame 7 Mm is no longer neglig i- 
ble (e.g., lHaas et al.ll200H iCharmandaris et al.l ll997L 
Spitzer IRS observations of bright radio-detected SMGs 
will be useful in quantifying the relationship between the 
5 — 8.5 /im and bolometric luminosities of these ultralu- 
minous sources. 

A relevant line of investigation is to determine what 
the various optical and near-IR color and magnitude se- 
lections imply for the infrared luminosity distributions 
of the galaxies they select. Figure shows the in- 
ferred Lir of galaxies with redshifts 1.5 < z < 2.6 
as a function of K s magnitude. We have assumed the 
IR/MIR flux ratio of 17.2 for the radio-detected SMGs 
shown in th e figure. Galaxies with LT S < 20 (e.g. , 
K20 samples: ICimatti et al.l2002al ICimatti et alJ2002bD 
have (Li r) ~ (1 — 2) x 10 12 L Q , similar to the value 
found bv lDaddi etUT po"05|) for K s < 20 BzK selected 
galaxies in the GOODS-N field. Alternatively, we find 
(Lir) - 5 x 10 11 L for U n GK, BzK, and DRG galaxies 
with 20 < K s < 21.0. As stated in § 15.11 the uncer- 
tainties on any individual value of log Lir are 0.3 dex for 
U n G7Z galaxies with spectroscopic redshifts and 0.6 dex 
for near-IR selected (BzK; DRG) galaxies with photo- 
metric redshifts. At any given K s magnitude, the range 
in Lir spans an order of magnitude assuming our con- 
stant conversion and larger than an order o f magnitude 
assum ing the two power-law conversion of lElbaz et alJ 
(2002). Finally we note that galaxies with K s < 20 at 
z ~ 2 which show some signature of star formation (i.e., 
those that are MIPS detected) generally have infrared 
luminosities that are a factor of two larger than those 
of galaxies with 20 < K s < 21. As discussed elsewhere, 
there is also a population of m assive galaxies with little 
detectable star formation (e.g.. Ivan Dokkum et alll200l 
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Fig. 8. — Le/i panel: Distribution of Lia as a function of _ft" s for galaxies in the various samples. We have assumed a IR/MIR flux ratio of 
17.2 for all galaxies including the radio-detected SMGs. Right panel: Lm as a function of (z — A')ab color. Symbols are the same as in the 
left panel. The arbitrarily normalized histogram indicates the distribution in (z — -K")ab color of DRGs and BzK/PE galaxies undetected 
at 24 fim. 



IReddv et al.ll200^ . 

We investigate this quiescent population of massive 
galaxies in the context of their star-forming counter- 
parts by examining Ltr as a function of (z — K)ab color 
(Figure (SJd) . The (z — K)ab color probes the Balmer 
and 4000 A breaks for galaxies at the redshifts of in- 
terest here, and i s also sensitive to the current star for - 
mation rate (e.g. JReddv"eT^l2005trDaddi et al.ll200l . 
Figure IHb shows that galaxies with redder (z — K)ab 
color (up to (z — K)ab ~ 3) have higher inferred Lir 
(and larger SFRs if the bolometric luminosity is at- 
tributed to star formation) on average than galaxies with 
bluer (z — K)ab colors, a trend similar to t hat found 
when stacking X-ray data ijReddv et alJl2005|l . A more 
interesting result is indicated by the histogram which 
shows the distribution in (z — K)ab color of DRGs and 
BzK/PE galaxies that are undetected at 24 /Ltm. Of 
the 13 MlPS-undetect ed DRGs and BzK/PE galaxies, 
7 have (z - K) AB > 3. IReddv et all lj2005|) found DRGs 
with (z — K)ab ^ 3 to have little X-ray emission and had 
colors si milar to those of IRAC Extremely Red Objects 
QEROs: lYan et al.ll2005j) . The lack of 24 /jm detections 
for these red (z — if )ab sources further supports the no- 
tion that they have little current star formation. It also 
rules out the possibility that they harbor Compton-thick 
obscured AGN as an explanation for their lack of X-ray 
emission, since we would then expect them to be bright 
at 24 /im. 

5.3. Stacked 24 \im Flux 

The high quality and depth of the MIPS data enable us 
to study the average properties of galaxies that are (in- 
dividually) undetected at 24 /im by stacking their emis- 
sion. The data were stacked by considering all galax- 
ies undetected at 24 ^m and X-ray wavelengths. To 
avoid contaminating the stacked signal, we only added 
galaxies to the stack if there were no bright MIPS 
sources within ks 5" of those galaxies. To ensure a 
clean signal, we extracted sub-images around each un- 




FlG. 9.— Stacked 24 /j,m flux of 48 isolated MlPS-undetected 
U n GlZ galaxies (/24/jm ^ 8 fiJy) with spectroscopic redshifts 1.5 < 
2 < 2.6, indicating a mean flux per undetected galaxy of (/24uro) ~ 
3.30 ± 0.48 nJy. 



detected galaxy, subtracted all detected point sources 
within those sub-images using the empirical PSF, and 
added the sub-images together. The stacked flux was 
measured in a manner similar to the measurement of 
detected MIPS sources. Figure shows the stacked im- 
age of 48 MlPS-undetected U n GTZ galaxies with spec- 
troscopic redshifts 1.5 < z < 2.6. The mean flux per 
object is (/24/im) ~ 3.30 ± 0.48 ^Jy, where the error 
is the dispersion of the background in the stacked im- 
age. At a mean redshift of (z) = 2.05, this flux cor- 
responds to Ltr w 2 x 10 10 Li©. Combining this re- 
sult with the mean L\n of MlPS-detected U n GlZ galax- 
ies implies a mean across the entire sample, neglecting 
"confused" sources, of w 2.3 x 10 11 L Q . This mean 
value does not change significantly if we add in directly- 
detected X-ray sources (including radio-detected SMGs 
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to Sg50/iin ~ 5 rnJy) because of their small number com- 
pared to the typical (less luminous) U n GlZ galaxy. The 
mean value of Lir ~ 2.3 x 10 11 L© is in excellent agree- 
ment with the average of Lir w 3 x 10 11 L Q found 
from st acked X-ray /radio and dust corrected UV esti- 
mates l|Reddv fc Ste idell 12004^ . This suggests that the 
non-detection of galaxies at 24 /im is due to them having 
lower SFRs and not because they are deficient in PAH 
emission for a given Lir. The advantage of the 24 /im 
data over X-ray/radio data is that we can estimate bolo- 
metric luminosities for individual L* (LIRG) galaxies at 
z ~ 2 rather than ensembles of galaxies. 

Combining our estimate of the MlPS-inferred average 
IR lumin osity of U n GTZ galaxies with the stacked radio 
results of lReddy fc Steidel |.2004), we find that the radio- 
IR relation appears valid on average for the sample. To 
quantify the radio- IR ratio for the sam ple, we compute 
the "q" parameter l)Condon et alJll99l|) : 



log( 



FIR 



3.75 x 10 12 W m- 2 ' 



log( 



Si 



4 GHz 



W m~ 2 Hz 



where Si. 4 GHz is the rest-frame 1.4 GHz flux density in 
units of W m -2 Hz -1 and 

FIR ee 1.26 x 10- 14 (2.58S 60 Mm + 5 100 „m) W m - 2 ,(4) 

where Seo ,, m and Si on ».m are the IRA S 60 and 100 /im 
flux densities in Jy ( He lou et ai1ll988[) . The implied "q" 
value for the U n GTZ sample is (q) ~ 2.5 if we assume 
log [560 pm/Sioo fj,m] ~ 0.2. This value of q is in excellent 
agreement with the value of q ~ 2. 4 found for > 10 11 L Q 
galaxies in the IRAS 2 Jy sample i)Yun et al.ll2001() . We 
also stacked the 24 /im data for undetected BzK/PE and 
DRG galaxies in the same manner as described above, 
which yielded a mean flu x pe r object of (/24/im) ~ 2.72 ± 
1.65 /iJy. As noted in § 15.21 most of these sources have 
very red (z-K)ab colors, and their low-level mid-IR and 
X-ray emission indicate they have low SFRs. Galaxies 
with /24/im ^$ 8 /iJy are discussed further below. 

6. DUST ATTENUATION IN OPTICAL AND NEAR-IR 
SELECTED GALAXIES 

Aside from inferring the infrared luminosity distribu- 
tions, we can use the MIPS data to examine the re- 
lationship betwe en dust extin c tion a nd rest-frame UV 
spectral slope. iMeurer et alJ Jl999) found a relation 
between the rest-frame UV spectral slope, j3, and the 
attenuation of local UV-selected starburst galaxies, pa- 
rameterized by the ratio Lfir/ L1600, where £1600 is the 
rest-frame 1600 A luminosity uncorrected for extinc- 
tion. This relation appears to fail, howe ver, for the most 
lumin ous starbursts such as ULIRGs (Goldade^^TJ 
12002ft and radio-detected SMGs {Chapman et alJ l2005l. 
A greater proportion of the star formation in galaxies 
with Lin ^ 10 12 L Q will be obscured by dust as com- 
pared with LIRG-type starbursts (e.g., see §01. There- 
fore, whatever UV emission is able to escape from the 
optically-thin regions of ULIRGs will constitute a lower 
fraction of the total bolometric luminosity. As a result, 
the rest-frame UV light can substantially underpredict 
(by a factor of 10—100) the bolometric luminosities of the 
most luminous starbu rsts, such as radio-detected SMGs 
ijChanman et alJl2005fl . Normal ( "quiescen t" ) star form- 
ing ga laxies also appear to deviate from the lMeurer et al.l 
( 1999) relation, such that Lfir/Li6oo is lower for a given 



amount of UV reddening tha n in starburst galaxies (e.g, 
lLaird et alJ 120051: iBelll 120021) a result that may be tied 
to the varying ratio of current to past-aver age star for- 
matio n rate of normal star forming galaxies ijKong et alJ 
2004). Alternatively, the star formation in local quies- 
cent galaxies (i.e., those with low SFRs) is more dis- 
tributed than in local starbursts so that a failure of the 
starburst reddening law may reflect a different distribu- 
tion of dust with respect to the star forming regions in 
quiescent galaxies. Observations of radio-detected SMGs 
and quiescent l y star -form ing galaxies suggests that the 
IMeurer et all l|1999ft and iCalzetti et all ( 2000) laws do 
not apply to these sources. 

The rest-frame UV spectral and mid-IR data of z ~ 2 
galaxies allow us to investigate how well the high redshift 
galaxies follow the local dust attenuation relation. The 
full SED modeling of U n GTZ galaxies in the GOODS-N 
field yields estimates of the best-fit star formation history 
(r), ag e, mass, SFR, and E(B — V) color exc ess for each 
galaxy l|Erb et all2006bllShar)lev et al.l2005D . The mean 
fractional uncertainty in E(B — V) is (awB-v)/E(B — 
V)} = 0.7 as determined from Monte Carlo simulations. 
To convert E(B — V) to /3 we assumed that 1 mag of 
extinction at 16 00 A (^4i6nn = 1) co rresponds to E(B — 
V) w 0.092 (e.g.. ICalzetti et al.120 001. For most galaxies, 
the best-fit star formation history is close to that of a 
constant star formation history (with decay time-scale 
r = 00). The most massive galaxies at z ~ 2 (with 
stellar masses M* > 10 11 M Q ) are generally better fit 
with declining star formation histories. We have assumed 
a CSF model for galaxies unless such a model provides 
a much poorer fit to the observed data than a declining 
star formation history. 

6.1. Results for Optically Selected Galaxies 

Figure EH shows £fir/-^1600 versus (5 for 
spectroscopically-confirmed U D G1Z galaxies with 
rcdshifts 1.5 < z < 2.6. The FIR luminosity is computed 
from £5-8.5^111 using the procedure described in § 15.11 
We estimated the rest-frame 1600 A luminosity from 
either the U a , G, or 1Z magnitude depending on the 
redshift of the galaxy. The majority of U n GlZ galaxies 
with inferred a ges > 100 Myr (soli d circles) appear to 
agree well with IMeurer et aLl l)1999j) relation, shown by 
the solid curve 12 . 

U n GlZ galaxies with the youngest inferred ages (< 
100 Myr; open circ les) exhibit a large offset from the 
IMeurer et alJ l)1999t) relation in the sense that they ex- 
hibit redder UV colors for a given dust obscuration than 
older galaxies which do follow the relation. Note that 
we have assumed a CSF model for the young galaxies 
shown in Figure ^| The inferred ages of these galax- 
ies are typically smaller than 50 Myr, which is approxi- 
mately the dynamical time across the galaxy. Assuming 
such small (and unrealistic ages) will cause us to over- 
estimate E(B — V) for these sources. The change in 
E(B — V) that results from fixing the age of the young 
galaxies to 100 Myr (AE(B - V) = 0.09) is not enough 
to complete ly account for the of fset of the young galax- 
ies from the IMeurer et al.l l|1999j) relation. This suggests 

12 Assuming the two power-law conversion to c ompute Lpm, 
result s in a similar distribution of points around the Meu rer et alJ 
(1999) relation but with larger scatter. 
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Fig. 10. — {Left) Dust absorption, parameterized by -FFiR./-/*igoo, versus rest-frame UV spectral slope, /3, for galaxies with redshifts 
1.5 < z < 2.6. Filled and open symbols respectively denote galaxies with inferred ages of > 10 Myr and < 100 Myr. The large pentagon 
shows the results for U n GlZ galaxies undetected at 24 /im, using the stacked results of § 15.31 The horizontal and vertical error bars on 
this point reflect the dispersion about the mean /3 and mean Lfir/Li600 of 24 fim undetec ted sources. The gr een points represent the 
results from an X-ray stacking analysis as described in the text. The solid line indicates the lMeurer et all I199 S) relation found for local 
UV-selected starburst galaxies. (Right) Same as left panel, but excluding galaxies younger than 100 Myr and color-coding objects by Lbol- 



that iCalzetti et al.l l)2000|) law may not be applicable to 
these young galaxies because of a different configuration 
of dust with respect to the star-forming regions, in which 
case a grayer extinction law may be appropriate. As one 
example, the well-studied lensed Lyman Break Galaxy 
MS1512-cB58 has an inferred age of - 70 - 100 Myr and 
millimeter continuum observations suggest that its in- 
frared luminosity i s smaller than one would predict from 
its UV reddening (Ba ker et al.l l2001). Regardless of the 
assumed extinction law, these young galaxies in the sam- 
ples examined here have similar bolometric (sum of ob- 
served IR and UV) luminosities as older galaxies (see 
§0- 

The deep X-ray data in the GOODS-N field allow us 
to estimate (X-ray inferred) average infrare d luminosities 



for well-defined samples of galaxies (e.g. , |Reddv_et_al 
| 2005tlLaird et al.l2005HNandra et all2002t Brandt et al 
2001). The green points in Figure UHl show the expected 
dust absorption inferred from the X-ray data as a func- 
tion of (3. These points were determined by stacking the 
X-ray data for non X-ray detected U n G7Z galaxies (with 
ages > 100 Myr) in bins of f3. We only considered stack- 
ing galaxies with ages > 100 Myr since the X-ray emis- 
sion is sensitive to the star formation rate once O and B 
stars evolve to produce high mass X-ray binaries, which is 
roughly 10 s years after the onset of star formation. The 
average X-ray flux per bin w as converted to a FIR flux 
using the lRanalli et alJ l|2003fl relation. Dividing the av- 
erage FIR flux per bin by the average 1600 A luminosity 
of objects in each bin yields the green points. The hori- 
zontal error bars indicate the bin width and the vertical 
error bars indicate the dispersion in stack ed X-ray flux 
(comp uted using the procedure des cribed in lReddv et al.l 
( 2005) and lReddv fc Steidefl l|2004ln added in quadrature 
with the dispersion in th e FIR/X-ray luminosity relation 
found bv iRanalli et all (120031) The X-ray results re- 
produce the iMeurer et al.l l)1999|) relation very well, pro- 
viding an independent confirmation that typical U n GlZ 



galaxies abide by the local dust obscuration relation for 
starburst galaxies. The X-ray and mid-IR data indicate 
that the UV slope can be used to deduce the extinction 
corrections for these typical galaxies and that such galax- 
ies have UV light that is only moderately extinguished in 
most cases. The agreement between the MIPS inferred 
FIR luminosities and X-ray inferred FI R luminosities 
(obtai ned with the empirically-d etermined |R analli et alJ 
( 2003) relation) suggests that the lRanalli et all l)2003|) ~ 
lation provides a reasonable description for most of the 
galaxies con s idered here. For comparison, adopting the 
iPersic et alJ l|2004|) relation would yield FIR luminosities 
five times in excess of those predicted from the 24 /im 
fluxes for the vast majority of z ~ 2 galaxies. 

We compute the average dust obscuration of U n GlZ 
galaxies undetected at 24 /im using the stacking results 
of § 15.31 and the result is denoted by the large pentagon 
in FigureUHl The stacked 24 fim FIR luminosity of MIPS 
undetected galaxies is divided by their average rest-frame 
1600 A luminosity. For these galaxies, £i6oo is on aver- 
age 1.6 times larger than their inferred Lfir indicating 
they are significantly less obscured than galaxies detected 
at 24 /im. These undetected galaxies also have relatively 
blue rest-frame UV SEDs (as indicated by their aver- 
age f}) compared to 24 /im detected galaxies. Further- 
more, the results of Figure ^indicate these faint sources 
have lower X-ray emission than 24 /im detected galaxies. 
All of these observations combined suggest that galaxies 
are undetected at 24 /im because they have lower SFRs 
and are less obscured than galaxies with brighter 24 /im 
fluxes. We further explore the nature of these MIPS un- 
detected sources in § [HI 

6.2. Results for Near-IR and Submillimeter Selected 
Galaxies 

Also shown in Figure ITU1 are BzK /SF galaxies to K s = 
21 not satisfying the U n GlZ criteria. As pointed out in 
S 15.21 these BzK /SF selected galaxies have inferred Lir 
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that are comparable to those of U a GlZ galaxies to the 
same K a limit, but of course with redder G — 7Z colors 
and a TZ ~ 0.5 mag fainter on average than U n GlZ galax- 
ies to the same K s limit. The results of FigurelTUIsuggest 
that BzK /SF galaxies lying outside the c olor range se- 
lected by the U n GlZ criteria also follow the lMeurer et al.l 
( 1999) relation. Similar to the results found in § 16. II 
for most U n G1Z galaxies, the mid-IR data indicate that 
the UV light from most BzK /SF galaxies appears to be 
moderately extinguished and that the UV slope can be 
used to estimate their attenuation. 

Almost all of the radio-detected SMGs considered 
here have inferred dust absorption factors (when we as- 
sume the 850 /im-inferred bolomctric luminosities) that 
are at least a magn itude larger than predicted by the 
iMeurer et alJ l)1999fl relation for a given rest-frame UV 
slope. The discrepancy is not as substantial (i.e., it is 
reduced by a factor of 2 — 10) if we inferred Ltr of the 
SMGs from their 24 /im fluxes assuming our conversion 
between MIR and IR luminosity. The U n G7Z criteria are 
designed to select objects where followup spectroscopy is 
feasible, and this usually implies setting a limit to the al- 
lowed E(B — V) (or (3) of objects in the sample. However, 
given that at least half the galaxies with Lboi ^ 10 12 L© 
have UV slopes comparable to that of the typical U n GTZ 
galaxy, it is not uncommon to find such dust-obscured 
galaxies in optical surveys. 

Of the limited sample of DRGs with photometric red- 
shifts 1.5 < z < 2.6, at least half lie above the lo- 
cal starburst attenuation l aw. We are able t o find 
DRGs that agree with the IMeurer et a l. (1999) rela- 
tion since the MIPS data studied here are significantly 
de eper (by a facto r of ~ 5) than the data considered 
in iPanovich et alJ l)2005t) . In particular, we find the 
surface density of DRGs between 1.5 < z < 2.6 with 
1 < log(F F i R /Fi 6 oo) < 2 of > 0.14 arcmin- 2 (this is 
a lower limit since there are number of DRGs without 
photometric redshifts, some of which may truly lie at 
redshifts 1.5 < z < 2.6), which is at least a factor of 20 
higher than in IPapovich et al.l l|2005|) . Our results sug- 
gest that the DRG population consists of galaxies with 
a very wide range in star formation rate, from galax- 
ies with little or no star formation (DRGs with very red 
(z — K)ab colors; § I5.2JI to those which are heavily ob- 
scured and rapidly forming stars. 

6.3. Relationship between (3 and Obscuration as a 
Function of Luminosity 

Figure ITUb shows galaxies with ages > 100 Myr within 
the samples, color-coded by their Lboi- Virtually all ob- 
jects with L bo i in the range 10 11 < L bo i < 10 123 L have 
j3 which appear to reproduce t heir obscuration as inferred 
from the lMeurer et all l)1999fl and lCalzetti et all l|2000ft 
laws. There is some weaker evidence that the galaxies 
with the lowest SFRs (undetected at 24 /im) as indicated 
by the green pentagon in Figure IlUb follow a different 
extinction law. More pronounced, however, is the sys- 
tematic offset of the most luminous galaxies considered 
here with Lboi > 10 12 3 L©. These ultraluminous galax- 
ies have rest-frame UV slopes that underpredict their 
obscuration by a factor of 10 — 100. The main results 
of Figure ITUb indicate that the relationship between UV 
reddening and obscuration is strongly dependent on the 
bolometric luminosity, but that most LIRG galaxies at 



z ~ 2 follow the local relation. 

7. RELATIONSHIP BETWEEN DUST OBSCURATION AND 
BOLOMETRIC LUMINOSITY 

The bolometric luminosity of star-forming galaxies can 
be well-approximated by the sum of the IR and UV lu- 
minosities as indicated in Equation El Figure ITT1 shows 
Lboi as a function of dust obscuration for objects in the 
various samples assuming a constant conversion relation 
between mid-IR and total IR luminosity. Typical (LIRG) 
galaxies at z ~ 2 will have Lboi ~ -^ir where ~ 90% of 
the bolometric luminosity is emitted in the infrared. The 
bolometric luminosity is strongly correlated with dust 
obscuration: galaxies with larger bolometric luminosi- 
ties are more dust obscured than less luminous galaxies. 
The best-fit linear trend for spectroscopically-confirmed 
U n GlZ galaxies detected at 24 /im is 

log Lboi = (0.62 ± 0.06) log -p^- + (10.95 ± 0.07) (5) 

Li6oo 

(solid line in Figure m we note that the two axes are 
not independent of each and may partly account for the 
tight scatter in the correlation). U n GlZ galaxies unde- 
tected at 24 /im are indicated by the pentagon. These 
undetected galaxies have an average bolometric luminos- 
ity of (Lboi) ~ 6 x 10 10 Lq and have UV luminosities 
that are a factor of ~ 10 less dust obscured than the typ- 
ical 24 /jm detected U n GlZ galaxy. Approximately half 
of the bolometric luminosity of these 24 /im undetected 
galaxies is emitted in the UV. Galaxies with inferred ages 
< 100 Myr (yellow symbols in Figure ITT|) have Lboi com- 
parable to those of older galaxies at z ~ 2, suggesting 
that these young galaxies have similar IR/MIR ratios 
as older galaxies. Therefore, the devia tion of the young 
galaxies from the IMeurer et al.l l)1999|) law as noted in 
§ suggests that we have over-estimated E(B — V) for 
these young sources and/or they may follow a steeper 
(e.g., SMC-like) extinction law. 

We note that U n G7Z galaxies with K s < 21 have bolo- 
metric luminosities and dust obscuration comparable to 
BzK galaxies to K s = 21 that are not optically selected. 
This implies that objects missed by optical selection but 
which appear in the near-IR selected BzK sample are 
not missed because they are more dust obsc ured, a result 
corro borated by X-ray stacking analyses ijReddv et alJ 
2005). The bright radio-detected SMGs have the highest 
inferred bolometric luminosities among all galaxies con- 
sidered here, of order ~ 10 13 L Q with dust obscuration 
factors > 100. Galaxies in the U n GlZ, BzK, and radio- 
detected SMG samples detected at 24 /im mostly follow 
the linear relation denoted by the solid line in Figure ITT1 
For SMGs, Lboi ~ Lir, so assuming the submillimeter es- 
timates of Lir (rather than the mid-IR estimates shown 
in Figure ITTJl will move the SMGs in a direction parallel 
to the z ~ 2 trend. 

DRGs detected at 24 /im span a large range in 
About half the DRGs follow the linear trend 
established for optically selected galaxies at z ~ 2. The 
remaining half of DRGs have similar bolometric lumi- 
nosities to those which follow the z ~ 2 trend, but the 
UV luminosities are a magnitude more at tenuated than 
what we would have predicted from the IMeurer et alJ 
(1999) law. The SED analysis (§ EJl demonstrates that 
all of the DRGs which follow the z ~ 2 trend are all rela- 
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Fig. 11. — Bolometric luminosity, approximated as the sum of the IR and UV luminosities, versus dust obscuration. Shown are the 
distributions for 2 ~ 2 galaxies, and the solid line indicates the best-fit linear relation for spectroscopically confirmed U n GlZ galaxies 
detected at 24 /ivci. For comparison, the pentagon denotes the result for galaxies in the U U G1Z sample und etected at 2 4 /im, and the yellow 
symbols denote objects with inferred ages < 100 Myr. We also show local galaxies from the compilation of Bell (2003). Galaxies of a given 
bolometric luminosity are typically 8 — 10 times less dust obscured at high redshift than at the present epoch. 



tively young galaxies (ages < 2 Gyr) and have lower stel- 
lar masses (M* < 10 11 M Q ). In contrast, DRGs which 
are offset from the trend are all older (ages > 2 Gyr) and 
all have masses > 1.2 x 10 11 Mq. The offset could be ex- 
plained naturally if the dust masses of galaxies increase 
as they age, a natural consequence of star formation. 
Note that if the massive, metal-rich DRGs have stronger 
PAH flux for a given IR luminosity than the younger 
galaxies, then this would serve to only increase the offset 
between the massive DRGs and the z ~ 2 trend. In fact, 
stacking the X-ray data for the younger and older DRGs 
indicates they have very similar bolometric luminosities, 
confirming the results obtained by inferring Lir from the 
PAH flux. The results of Figure ^2 suggest that much 
of the dust in galaxies with the largest stellar masses 
was produced by star formation prior to the episode cur- 
rently heating the dust. Therefore, such galaxies will 
have larger dust obscuration for a given bolometric lu- 
minosity. Assuming the two power-law conversion of 



lElbaz et all {2002) would result in a z ~ 2 trend with 
a slope 20% larger than given in Equation |5j but with 
approximately the same intercept within the uncertain- 
ties, so our conclusions would be unchanged. 

As galaxies are enriched with dust as they age, then we 
expect to see an even greater difference in dust obscura- 
tion between z ~ 2 galaxies and those at the present 
epoch. To investigate this, we examined Lboi versus 
Ltr/ . Li piop for the sample of local galaxies compiled by 
iBelll (|2003f) . shown by crosses in Figure 11 ^ The local 
sampl e includes the ULIRGs studied by Gol dader et alJ 
(2002). Unfortunately, the UV and IR data for local 
LIRG and ULIRGs are relatively sparse. However, of 
the small sample of local galaxies with Lboi ^ 10 11 L©, 
almost all (10 of 11) lie to the right of the linear trend 
at z ~ 2 and at least half occupy the same region as the 
old, massive DRGs at z ~ 2. In fact, an interesting corol- 
lary to the above discussion is that massive, star-forming 
DRGs at z ~ 2 are more analogous to local ULIRGs than 
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bright SMGs at z ~ 2, both in terms of bolometric lu- 
minosity and dust-obscuration. Local ULIRGs undoubt- 
edly carry a significant amoun t of dust into their cur rent 
star formation episodes (e.g., iGoldader et alJ 12002)) . so 
it not surprising that they have similar dust-obscuration 
factors as massive, star-forming DRGs at z ~ 2. On the 
other hand, most galaxies lying on the z ~ 2 trend, in- 
cluding many bright SMGs, are likely undergoing their 
first major episode of star formation and have relatively 
low dust-to-gas ratios, unlike the more massive (offset) 
DRGs and local galaxies. 

The offset between z = and z ~ 2 galaxies can be 
seen at fainter bolometric luminosities where the local 
sample includes more galaxies (10 10 < Lhn i 5, 10 11 ). Re- 
stricting our analysis to galaxies in the Bell ( 2003) sample 
with Lboi comparable to those of 24 /an undetected z ~ 2 
galaxies, we find that the local sample is on average ~ 10 
times more dust obscured than 24 /im undetected galax- 
ies at z ~ 2. Further, recent GALEX results indicate 
that local near-UV selected galaxies with Lboi = 10 11 L 
have a mean dust obscuration factor of ~ 10; this is 
8 times larger than the inferre d dust obscuration of a 
Lboi = 10 11 L galaxy at z = 2 ijBurgarella et al.l f2005K 
To summarize, the important result from Figure ^JJ is 
that galaxies of a given bolometric luminosity are on av- 
erage a factor of 8 — 10 less dust obscured at z ~ 2 than 
at the present epoch, confirm ing the trend first noted by 
lAdelberger fe Steidell (|2000f) between galaxies at z = 0, 
z ~ 1, and z ~ 3; this result is also suggested by the work 
of lGalzetti fc Heckmanl l)1999l) . A gain, this result could 
be anticipated if successive generations of star formation 
add to already existing dust within galaxies and/or if the 
dust distribution within galaxies becomes more compact 
with time (e.g., via the effects of mergers which tend to 
drive gas and dust to the central kpc of the system). 
The net result of dust enrichment and a more compact 
distribution of dust (e.g., after a merger event) is an in- 
crease in the dust column density towards star-forming 
regions. The relationship between dust obscuration and 
Lboi (i-e., EqEJ indicates that for the mean Lboi of the 
U n GlZ selected sample of galaxies of Lboi ~ 2.3x 10 11 L , 
the average dust obscuration is (Lir/Li 6 oo) « 4.6. This 
factor is in excellent agreement with the mean attenu- 
ation of 4.5 — 5.0 infer red from stacked X-ray analyses 
ijReddv fc Steidell l2004). One would observe a factor of 
4 — 5 attenuation in a galaxy one order of magnitude less 
luminous at z — than at z ~ 2. The implication is 
that, while it is true that a larger fraction of star forma- 
tion at high redshifts occurs in dustier systems, the dust 
obscuration we observe for galaxies of a given Lboi has 
less of an impact on observations of high redshift galax- 
ies than one would have surmised on the basis of present 
day galaxies. 

8. PROPERTIES OF 24 FAINT GALAXIES 

8.1. Ages and Masses of Faint 24 [im Galaxies 

In addition to the information gleaned from the stack- 
ing analysis described above, we also have detailed infor- 
mation on the stellar populations of galaxies with faint 
mid-IR emission. Optical and near-IR selected galaxies 
undetected at 24 /im appear to have a distribution in ages 
which is similar to that of 24 /im detected galaxies (left 
panel of Figure I12|) , so young ages cannot explain why 
they are undetected at 24 /im. Alternatively, although 



we find a large range in inferred stellar mass of galaxies 
with /24/im < 8 /iJy, the mean stellar mass of undetected 
objects is 0.4 dex lower in logM* than 24 /mi detected 
galaxies (right panel of Figure [T^ . Regardless of these 
small differences in the stellar populations of 24 /im de- 
tected and undetected sources, the X-ray data indicate 
that the primary reason why galaxies are undetected at 
24 /im is because they have lower SFRs (Figure We 
demonstrate in the next section how the rest-frame UV 
spectral properties of galaxies can be used to interpret 
their 24 /im emission. 

8.2. Composite UV Spectra 

A unique advantage of our optical U n GlZ selection is 
the efficiency with which we are able to obtain rest-frame 
UV spectra for these galaxies, and this spectroscopy al- 
lows for an independent probe of the physical conditions 
in the ISM. While the S/N of any individual spectrum is 
typically too low to accurately measure interstellar ab- 
sorption line widths, we can take advantage of the large 
number of spectra by stacking them to create a higher 
S/N composite spectrum. To investigate differences in 
the ISM as a function of infrared luminosity, we con- 
structed composite UV spectra for (a) the top quartile 
of 24 /im detected U n GlZ galaxies, and (b) all U n GTZ 
galaxies undetected at 24 /im. In order to stack the spec- 
tra, we first de-redshifted them by the systemic redshift. 
The systemic redshift was inferred from a weighted com- 
bination of the measured absorption an d /or emission line 
redshi fts, following th e procedu r e of l Adelberger et all 
( 2003). We used the iRix et alJ 1)20041) prescription to 
normalize the composite spectra to the underlying stel- 
lar continua. The detected and undetected composite 
spectra consist of 39 and 73 galaxies, respectively, and 
are shown in Figure 1131 For comparison, the mean 
24 /tm flux of MIPS detected and undetected galaxies 
is (/24 fim) ~ 100 /iJy and ~ 3 /iJy, respectively; both 
sub-samples have (z) ~ 2.1. Table lists the measured 
rest-frame equivalent widths of several interstellar ab- 
sorption lines in the composite spectra. The primary 
difference between the rest-frame UV spectra of 24 /im 
detected and undetected galaxies is that the latter have 
interstellar absorption lines that are a factor of 2 weaker 
than the lines in the 24 /im detected galaxies. Because 
the line strengths are controlled by the combination of 
the velocity spread in outflowing interstellar material and 
the covering fraction of optically-thick material, this in- 
dicates that galaxies weak in mid-IR emission are likely 
to have more quiescent ISM than 24 /im-bright galaxies, 
a result expected if those galaxies undetected by MIPS 
have lower SFRs, and hence lower energy input into the 
ISM and a lower level of dust enrichment, than 24 /im 
detected galaxies. 

Comparing the mid-IR detections with nondetections, 
we find the latter have significantly stronger Lya emis- 
sion than the former. The emergent Lya profiles of galax- 
ies will depend strongly on a number of physical param- 
eters including the spectrum of UV radiation (i.e., the 
stellar IMF), presence of outflows, and dust covering frac- 
tion. Neglecting all of these effects, galaxies with larger 
SFRs will have stronger Lya emission. However, given 
that the bolometric luminosity of star-forming galaxies 
scales with dust-obscuration (e.g., Figure ITT)) . we might 
expect 24 /tm detected galaxies to have larger dust col- 
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Fig. 12. — Histograms of the age (left panel) and mass (right panel) distributions for 24 fim detected and non-detected U n GTZ galaxies 
with redshifts 1.5 < z < 2.6. The age distributions of detected and non-detected U n GlZ galaxies are similar. The mean of the mass 
distributions are offset such that undetected galaxies have (logM*) that is 0.4 dex lower than that of 24 (im detected galaxies. 
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Fig. 13. — Normalized composite UV spectra for the 39 galaxies in the top quartile of 24 (im detected U n G7Z galaxies (top) with 
(/24fi.m) ~ 100 fi-]y and 73 !/„GR galaxies undetected at 24 /an (bottom) with (/24um) ~ 3 /^Jy. 
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TABLE 3 

Interstellar Absorption Line Wavelengths and 
Equivalent Widths for 24 /jm Detected and 
Undetected U n G1Z Galaxies 



Line 


\ a 
^trans 


Detected W A b 


Undetected W A b 


Si II 


1260.4 


2.58 


0.99 


OI + Si II 


1303.3 


2.27 


1.58 


C II 


1334.5 


2.05 


0.94 


Si II 


1526.7 


2.36 


1.32 


Fe II 


1608.5 


0.87 


1.15 


Al II 


1670.8 


1.85 


1.15 


Al III 


1854.7 


1.66 


0.55 


Al III 


1862.8 


0.81 


0.64 



a Transition wavelength in A. 

b Measured rest-frame equivalent width in A. 



umn densities than undetected galaxies; this may partly 
explain the absence of Lya emission in 24 /jm detected 
galaxies. In addition, the velocity spread of the ISM will 
also affect the emergent Lya profile: the larger veloc- 
ity spread in 24 /im detected galaxies, as indicated by 
their stronger interstellar absorption lines, implies most 
Lya photons will have larger scattering path-lengths and 
are more likely to be attenu ated by dust and /or scat- 
tered out of resonance (e.g.. lHansen fe Peng O hi 120051: 
lAdelberger et~aT]l200^ . 

Finally, we note that the stacked X-ray analysis of 
24 /zm-undetected galaxies confirms they have lower 
SFRs than 24 /xm-detected galaxies. Therefore, galax- 
ies are undetected at rest-frame 5 — 8.5 /xm primarily 
because they have lower SFRs and not because they are 
deficient in mid-IR PAH luminosity for a given Lfr. If 
such undetected galaxies had depressed MIR/IR flux ra- 
tios, we would not have expected to see as large a differ- 
ence in the strengths of their interstellar absorption lines 
as compared with 24 /xm detected galaxies. 

9. MID-IR PROPERTIES OF MASSIVE GALAXIES AT Z ~ 2 

The epoch between z = 3 and z = 1 appears to be 
the most active in te r ms of the buildup of stellar mass 
Ce.g.. lDickinson et al.1 1)2003|) . see also §P), but significant 
numbers of massive galaxies (M* > 10 11 L Q ) already ap- 
pear to be in place by rcdshifts z ~ 2. The subsequent 
evolution of these massive galaxies and their relation to 
the local population of massive and passively evolving 
elliptical galaxies is an important question. It is useful 
to determine, therefore, what the mid-IR properties of 
massive galaxies tell us about their bolometric luminosi- 
ties. 

Figure PHI shows the bolometric luminosity of galaxies 
in the U n GlZ, BzK, and DRG samples as a function of 
best-fit stellar mass. The mass modeling comes from the 
SED analysis where we have fit for the R+ J K S +IRAC 
photometry assuming a constant star formation (CSF) 
model. As discussed elsewhere, the stellar mass is the 
most well determined parameter from the SED analysis 
and is relati vely insensitive to the assumed star forma - 
tion history l)Shat)lev et al.ll2005l IPapovich et al.ll2001^ . 
There are several interesting aspects worth considering 
in Figure El First, we note that galaxies with the 
youngest inferred ages (< 100 Myr) have low stellar 
masses (M* < 2 x 10 10 M©) and span a large range in 
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Fig. 14. — Inferred Lt, i as a function of best-fit stellar mass 
assuming a CSF model for galaxies in the UnGlZ, BzK, and DRG 
samples. We assumed a IR/MIR ratio of 17.2 to convert L5_8.5/jm 
to Lir. The large solid pentagon and square denote results for un- 
detected U n G1Z and the 7 undetected DRG galaxies, respectively. 
The light gray (yellow) symbols indicate 24 ^tm-detected galaxies 
with inferred ages < 100 Myr. 

bolometric luminosity from LIRG- to ULIRG-type galax- 
ies. Among U n GlZ galaxies with spectroscopic redshifts 
1.5 < z < 2.6, we note a trend of increasing Lboi with 
stellar mass; U n G1l galaxies with M* > 10 11 M have 
infrared luminosities typical of ULIRGs. The BzK /SF 
and DRG criteria cull M* > 10 11 M galaxies with a 
large range in bolometric luminosity, from 10 11 L to 
> 10 12 L . Optically selected sources undetected at 
24 /xm have stellar masses similar to those of detected 
U n GlZ galaxies, but with Lboi a factor of ~ 10 lower. As 
noted before, the 24 /xm emission from these galaxies ap- 
pears to be primarily dependent on their current star for- 
mation rate. Finally, we note that DRGs undetected at 
24 /xm have similar stellar masses as those which are de- 
tected (M* > 10 11 M Q ), but are on av erage 40 times less 
bolometrically luminous. As noted in § 15.21 these 24 /im- 
undetected DRGs have very red (z — K)ab colors. Their 
faint 24 /xm emission corroborates the results of X-ray 
stacking analyses that indic ate these galaxies h ave very 
little current star formation l)Reddv et alJl200'5l . 

The results of Figure E] suggest that galaxies with 
large stellar masses at z ~ 2 span a large range in bolo- 
metric luminosity, from galaxies with little star formation 
to ULIRG-type systems. The range is likely larger than 
what one would infer from Figure 1141 since we have ex- 
cluded directly-detected X-ray sources that may be heav- 
ily star- forming galaxies and/or host AGN. Figure ITU 
suggests that > 70% of massive galaxies (M* > 10 11 M ) 
with redshifts 1.5 < z < 2.6 in optical and near-IR sur- 
veys have Ltd > 3 x 10 11 L ( or SFRs > 20 yr" 1 , 
assuming the SFR calibration of lKennicuttlll"99 8a ) . Our 
results imply that most 1.5 < z < 2.6 galaxies with large 
stellar masses (M* > 10 11 M ) have levels of star for- 
mation that exceed those of LIRGs. 

10. DISCUSSION 

10.1. Selection of LIRGs and ULIRGs at z ~ 2 

In § 15.21 we showed that optical and near-IR selected 
samples to 1Z = 25.5, or K s = 22, host galaxies with 
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a wide range in infrared luminosity, from a few times 
10 10 L Q up to the most luminous objects at these red- 
shifts with Lir > 10 12 L Q . Typical galaxies in these 
samples have infrared luminosities in the range 10 11 < 
Lm < 10 12 Lq, similar in luminosity to, but with dust- 
obscuration a factor of ~ 10 lower than (see § 0), lo- 
cal LIRGs. One advantage of optical (i.e., rest- frame 
UV) selection of LIRGs and ULIRGs at high redshift 
is that it allows for the kind of efficient spectroscopic 
followup necessary to accurately interpret the nature of 
24 /jm emission from these galaxies. As shown in § 
the if-correction depends sensitively on redshift such 
that even for "good" photometric redshift estimates with 
Sz/(1 + z) ~ 0.1, the corresponding uncertainty in the 
rest-frame 5 — 8.5 /im luminosity increases by a factor of 
5. A particularly unique aspect of our study combining 
the optical sample with MIPS observations is that the 
spectroscopic database can be used to assess the physi- 
cal conditions of the ISM in galaxies as function of I/bob 
providing an additional method for probing the detailed 
nature of 24 /jm galaxies at z ~ 2 (§ I8.2|l . Aside from 
the constraints on the mid-IR luminosities possible with 
spectroscopic redshifts, precise positions of sources from 
higher spatial resolution and shorter wavelength data en- 
able the deblending of most z ~ 2 galaxies. The de- 
blending procedure made possible by optical, near-IR, 
and Spitzer IRAC observations enable accurate identifi- 
cation and photometry of faint galaxies well below the 
MIPS 24 fim confusion limit and will provide a more 
complete "census" of the LIRG population at redshift 
z ~ 2 than possible using MIPS observations alone. 

Further, selection by optical colors gives important in- 
formation on the unobscured component of the star for- 
mation in galaxies and complements well the information 
on the obscured component probed by the 24 fim obser- 
vations. Objects with lower star formation rates will have 
bolometric luminosities that are typically dominated by 
the observed UV emission and objects with larger star 
formation rates will have bolometric luminosities that 
are dominated by the observed IR emission. The tran- 
sition between the UV and IR dominated regimes (i.e., 
where Lir = £1600) at z ~ 2 occurs for galaxies with 
Lir w 10 11 L Q , or about 0.3L* (Figure ITU . A compar- 
ison with the z — sample of Bell (2003) shows that 
the bolometric luminosities begin to be dominated by IR 
emission (i.e., £ir/£i6oo > 1) for galaxies which are two 
orders of magnitude more luminous at z ~ 2 than at the 
present epoch. As discussed in § this is plausibly ex- 
plained as a result of higher dust-to-gas ratios in the local 
galaxies. More generally, galaxies of a given dust obscu- 
ration are anywhere from 2-100 times more luminous at 
z ~ 2 than locally (Figure lll|) , with the greatest dif- 
ference for galaxies with relatively low £ir/£i600 % 20. 
The implication of these observations is that while it is 
certainly true that a larger fraction of the star forma- 
tion at high redshifts occurs in dustier galaxies, selec- 
tion via rest-frame UV colors (and performing followup 
spectroscopy) is easier at high redshift than locally for 
galaxies at a given bolometric luminosity. Optical selec- 
tion is therefore arguably the most promising and spec- 
troscopically efficient method for selecting LIRGs (which 
undoubtedly account for a significant fraction of the star 
formation rate density and far-infrared background; e.g., 
lAdelherger Steidefaofl) at z ~ 2. 



As demonstrated in § 15.21 and in iReddv et al.l {2005) , 
ULIRGs and SMGs at these redshifts also often ap- 
pear in optical and near-IR selected samples; ~ 50% 
of the most luminous SMGs have enough unobscured 
star formation that they satisfy the z ~ 2 optical crite- 
ria l|Ohapman et al.ll2005lK Accounting for the obscured 
portion of star formation in these ultraluminous sources 
of course requires the longer wavelength data since rest- 
frame UV slopes underpredict their Ltr (§EJ- ln some 
sense, 24 /im observations are a more powerful method 
of estimating Lboi of these ultra-luminous sources since 
their detection significance at 24 fim is typically 50 — 100 
times larger and the beamsize is a factor of 9 smaller than 
at 850 fim. The ability to uniformly cover large areas of 
the sky with MIPS observations provides an additional 
advantage over current submillimeter surveys which suf- 
fer from areal and depth incompleteness, thus making it 
difficult to constrain the volumes probed. Using 24 fim 
observations to assess the global energetics of ultralumi- 
nous sources of course requires that we accurately cali- 
brate the i5-8.5/jm/iiR ratios for these objects. 

We have demonstrated that the typical galaxy in op- 
tical and near-IR samples of z ~ 2 galaxies has I/ir cor- 
responding to that of LIRGs. A related issue is whether 
most LIRGs at z ~ 2 can be selected by their rest-frame 
UV or optical colors. A direct comparison of the number 
counts of MIPS sources to the number of 24 fim detected 
U n GTZ and BzK galaxies to 8 fiJy (the GOODS-N MIPS 
24 fim 3 a sensitivity limit) is not possible since (a) we 
primarily relied on the -fT s -band data to deblend sources 
in the 24 /im imaging and (b) the redshift distribution 
of MIPS sources to 8 fiJy is not yet well established. 
Nonetheless, including both optical and near-IR selected 
LIRGs at z ~ 2 ensures that we must be reasonably 
"complete" for both optically-bright (and 24 /xm-faint) 
LIRGs in the optical sample and optically-faint (and 
24 /xm-bright) LIRGs in the near-IR selected BzK and 
DRG samples. Galaxies with LIRG luminosities will pre- 
dominantly have near-IR magnitudes bright enough to 
be considered in our analysis (i.e., K s < 21 according to 
the left panel of Figure|BJ). Objects not selected by these 
various criteria will likely cither fall at different redshifts, 
not have LIRG luminosities, and/or may be scattered out 
of th e color selection wi ndows due to photometric error 
(e.g.. IReddv et aLl l2005). As an example of one form of 
photometric scatter, in the course of the z ~ 3 Lyman 
Break Galaxy Survey, we relied on U n GlZ photometry 
based on images of the HDF- North field taken at the 
Palomar Hale 5 m Telescope (jSteidel et alJl2003|) . We 
subsequently imaged a larger portion of the GOODS- 
North field using the Keck I Telescope l)Steidel et alJ 
2004). Our photometric analysis indicates that of the 
BX/BM objects to 1Z = 25.5 identified in the Palomar 
imaging, about 76% were recovered as BX/BM objects 
in the Keck imaging. A small fraction of the remain- 
ing 24% were recovered using LBG selection. The level 
of scatter between different photometric realizations in 
other fields is also typically ~ 25% and is mostly due to 
the narrow photometric windows used to select BX/BM 
galaxies. Regardless of these photometric effects, it is 
highly unlikely that there exist large numbers of LIRGs 
at z ~ 2 with such different optical and near-IR proper- 
ties that they would be completely absent from all the 
samples considered here. Finally, our knowledge of the 
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exact positions of optical and near-IR selected galaxies 
from the higher spatial resolution .ft" s -band and IRAC 
data allows us to mitigate the effects of confusion (see 
§ 01 , so we should be reasonably complete for galaxies 
which are detected at 24 /zm to 8 /iJy but which might 
otherwise be confused with brighter sources. It is there- 
fore reasonable to conclude that the LIRG population at 
z ~ 2 is essentially the same population of galaxies that 
are selected in optical and near-IR samples. 

10.2. Mass Assembly at High Redshift 

We demonstrated that LIRGs and ULIRGs are present 
over the full range of stellar mass, from ~2x 10 9 M to 
5 x 10 11 M Q , for galaxies in the samples considered here 
(FigureEJ- To assess the significance of the current star 
formation in the buildup of stellar mass, we have com- 
puted the specific star formation rate, <f>, defined as the 
SFR per unit stellar mass. We show the observed <f> for 
galaxies in our sample as function of stellar mass in Fig- 
ure El The correlation between and M* could have 
been predicted from Figure PHI since the range of iboi is 
similar over the range of M* considered here. We also 
note that the correlation is accentuated since (a) </> is not 
independent of M* and (b) there are presumably galaxies 
with low cj) and low M* that would be missing from the 
optical and near-IR samples (irrespective of the MIPS de- 
tection limit). Furthermore, the upper envelope of points 
in Figure 1151 is defined by our cut to exclude luminous 
AGN based on the Chandra X-ray data. Nonetheless, we 
find that star-forming galaxies with large stellar masses 
M* > 10 11 M Q without an AGN signature have specific 
SFRs that are 1 to 2 orders of magnitude lower than 
those of young galaxies (yellow symbols), implying that 
the current star formation contributes more significantly 
to the buildup of stellar mass in low mass galaxies than 
high mass galaxies at z ~ 2. This change in specific 
SFR as a function of mass has be en observed at later 
epochs as well (e.g.. lBell et a l. 2005). The shifting of the 
relationship to lower specific SFRs at later epochs for 
galaxies with l arge stellar masses has been referred to as 
"downsizing" l|Gowie et, al.lll99fift . 

The most massive, star-forming DRGs at these red- 
shifts with large dust-to-gas ratios (see Figure lll|) also 
have the lowest <j) for MlPS-detected galaxies at z ~ 2. 
IPapovich et alJ l|2005j) demonstrate that the integrated 
specific SFR of these massive DRGs at z ~ 2 is ~ 2 
magnitudes larger than for galaxies with similar masses 
(M* > 10 11 Mq) at lower redshifts 0.3 < z < 0.5, based 
on a c omparison with the COMBO- 17 sample ( Be ll et al.l 
12005ft . This decrease in suggests that M* > 10 11 M 
galaxies have built up most of their stellar mass by z ~ 2. 
In fact, the mass-doubling time becomes increasingly 
large for these massive galaxies if we assume the case 
of an exponentially declining star formation model. In 
this case, <j) will evolve with time t as 



(l + /)exp(-t/r) 
t[1 - exp(-t/r)] : 



(6) 



where / is the ratio of the gas mass that is lost due to 
outflows to that formed in stars and r is the star for- 
mation decay timescale as defined in § 0] Given some 
initial SFR, a CSF model (i.e., with r = oo) will track 
straight through the points corresponding to galaxies 



with Al* > 10 11 M©, but such a model would pre- 
dict a <fi at z ~ 0.4 that is ~ 5 times large r than ob- 
serve d at z ~ 0.4 in the COMBO-17 survey ( Be ll et al] 
2005) . Therefore, a declining star- formation history may 
be more appropriate for describing the future evolution 
of ga laxies with large s tellar masses at z ~ 2. In partic- 
ular, lErb et aD j2()()(ia ■ find that a model that assumes 
a super-solar yield of metals (i.e., ratio of mass of metals 
ejected into the ISM to mass of metals locked in long- 
lived stars) of y = 1.5 Z Q and an outflow rate of 4 x SFR 
(/ = 4) appears to best-fit the observed metallicities of 
U n G7Z galaxies at z ~ 2 as a function of gas fraction. If 
all star-forming galaxies at z ~ 2 follow a similar evolu- 
tionary track as the M* > 10 11 M© galaxies (i.e., follow 
an exponentially declining star formation history with 
large outflow rate), then the scatter of galaxies with a 
given specific SFR simply reflects the range in the final 
stellar ma sses and dark ma tter halo masses (see also dis- 
cussion in lErb et alJl2006ah . 

We can directly relate the specific SFR cf> with the cold 
gas fraction, fj, = M gas /(Mg as + M*). If we assume that 
the SFR is proportional to the cold gas mas s, M Fasi to 
the 1.4 power according to the Schmidt law (Kennicutt 
Il998b[) . then 



4> oc 



M 1A 

gas 

M* ' 



It then follows that 



4> = c 



1- M ' 



(7) 



(8) 



where C is a constant that depends on the constant 
of proportionality between the SFR and gas mass sur- 
face densities in the Schmidt law and the total gas mass 
at virialization (i.e., when star formation commences). 
There is a one-to-one correspondence between the spe- 
cific SFR, cf>, and gas fraction, /i, such that galaxies with 
large specific SFRs will have a larger fracti on of cold 
gas tha n galaxies with small specific SFRs. lErb et all 
(2006b) demonstrate that /i decreases as a function of 
stellar mass for a large sample of U n GlZ galaxies with 
Ha spectroscopy, with a mean \i across the sample of 
(fi) ~ 0.5. The mean specific SFR for these galax- 
ies is ((f)) ~ 3 Gyr -1 . Using these mean values to es- 
timate C, we show the range of fi on the right-hand 
axis of Figure El Young galaxies with ages less than 
100 Myr in our sample also have the largest gas frac- 
tions (ju ~ 0.6 — 0.9) and largest specific SFRs (Fig- 
ure El compared with older galaxies at z ~ 2. These re- 
sults strongly suggest that the young galaxies have large 
reservoirs of cold gas and have just begun forming stars. 
Similarly, galaxies with the largest stellar masses have 
lower specific SFRs and therefore lower cold gas fractions 
(fi ~ 0.1 — 0.3) and are likely to cease star- formation 
in a relatively short time. The correlation between <fi 
and M* revealed by Figure El is actually expected from 
the trend between /x and M* inf e rred fro m the Ha spec- 
troscopic analysis of lErb et alJ (2006b). Furthermore, 
lErb et all l)2006aft find a tight trend between metallic- 
ity and gas fraction such that galaxies with lower /i are 
more metal-rich. These galaxies will therefore be more 
dust-obscured since metallicity is directly proportional 
to dust-to-gas ratio. Our analysis confirms our expecta- 
tion that the most dust-obscured objects at z ~ 2 (i.e., 
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Fig. 15. — Specific SFR, <j>, as a function of stellar mass for galaxies in the samples considered here. Symbols are the same as in Figure H^I 
(yellow symbols indicate galaxies with inferred ages < 100 Myr). The right-hand axis indicates the gas fraction fi associated with galaxies 
of a given <j>. 



the DRGs with the largest dust-to-gas ratios, Figure ITTJI 
also have low gas fractions as demonstrated in FigurclT51 
Figure H5l indicates that the trend in gas fraction versus 
stellar mass (or, metallicity versus gas fraction ) found for 
U n GlZ -se\ected galaxies with Ha spectroscopy lErb et all 
2006a) also applies to galaxies selected by their near-IR 
colors. 

Combining the optical and near-IR samples, we find 
a wide range in the evolutionary state of galaxies at 
z ~ 2, from those which are just starting to form stars 
to those which have already accumulated most of their 
stellar mass and are about to become passive or already 
are. Almost all DRGs with the reddest (z — K)ab color 
((z — K)ab > 3) are undetected at 24 /jm and in the deep 
Chandra X-r ay data (e.g., see rig ht panel of Figure[B]and 
Figure 16 in iReddv etafl 12005). The specific SFRs of 
these red DRGs are the lowest observed for the z ~ 2 
galaxies considered here; they have cold gas fractions of 
less than 5%, confirming that they have essentially shut- 
off star formation by z ~ 2. The star- forming DRGs 
with M* > 10 11 M© (which are dusty and metal-rich 
as inferred from their large dust-to-gas ratios; cf., Fig- 
ure lll[) are likely to reach this passively-evolving state 
by redshifts z < 1.5. A simple model which assumes a 



high outflow rate proportional to the SFR and exponen- 
tially declining star formation history (i.e., the model 
which appears to best-fit the obse rved metallicities of 
galaxies as function of gas fraction; lErb et aTll2006aj) is 
sufficient to explain the presence of galaxies with large 
stellar masses and little star formation by z ~ 2 with- 
out invoking AGN feedback, despite the fact that a large 
AGN fraction of ~ 2 5% is observed am ong galaxies with 
large stellar masses l|Reddv et aUl2005|) . 

11. CONCLUSIONS 

We use Spitzer MIPS data to examine the bolomet- 
ric luminosities and extinction properties of optical and 
near-IR selected galaxies at redshifts 1.5 < z < 2.6 in 
the GOODS-North field. At these redshifts, the mid-IR 
(5 — 8.5 /im) features associated with polycyclic aromatic 
hydrocarbon (PAH) emission, which are ubiquitous in lo- 
cal and z ~ 1 star-forming galaxies, are shifted into the 
MIPS 24 fim filter. Extensive multi-wavelength data in 
the GOODS-North field, including very deep Chandra X- 
ray data, allow us to test the validity of inferring the bolo- 
metric luminosities of high redshift galaxies from their 
rest- frame 5 — 8.5 (im emission. Galaxies at z ~ 2 are 
selected by their optical (U n GlZ) and near-IR (BzK and 
J—K s ) colors, and for comparison we also consider radio- 
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detec ted submillimeter galaxies (SMGs; iChapman et al.1 
2005) . The optically-selected sample is advantageous be- 
cause we have knowledge of the precise redshifts of > 300 
galaxies from spectroscopy, allowing for the most accu- 
rate estimates of the rest-frame mid-IR luminosities of 
z ~ 2 galaxies. In addition to our spectroscopic sample, 
we use deep optical, near-IR, and Spitzer IRAC data to 
derive photometric redshifts for BzK galaxies and those 
with red J — K s > 2.3 colors (Distant Red Galaxies; 
DRGs). The principle conclusions of this study are as 
follows: 

1. Using local templates to iif-correct the observed 
24 fim fluxes, we find that the rest-frame 5 — 8.5 /mi lu- 
minosity (£5— 8.5/im) of z ~ 2 galaxies correlates well with 
their stacked X-ray emission. A subset of galaxies with 
Ha measurements have Ha-inferred bolometric luminosi- 
ties that correlate very well with their Ls-s.s^m-inferred 
bolometric luminosities (with a scatter of 0.2 dex). These 
observations suggest that L 5 - 8 . 5fJiln provides a reliable es- 
timate of Lir for most star forming galaxies at z ~ 2. 

2. We find that the optical and near-IR selected z ~ 2 
galaxies span a very wide range in infrared luminosity 
from LIRG to ULIRG objects. We find a mean infrared 
luminosity of (£ir) ~ 2 x 10 11 L Q for galaxies in the 
optical and near-IR samples, in excellent agreement with 
the value obtained from a stacked X-ray analysis. The 
optical and near-IR selected star-forming galaxies likely 
account for a significant fraction of the LIRG population 
at z ~ 2. Galaxies with K s < 20 have Ljr greater than 2 
times that of galaxies with K s > 20.5. Non-AGN galaxies 
with the reddest near-IR colors ((z — A")ab > 3) are 
mostly undetected at 24 //m suggesting they have low 
SFRs, a conclusion supported by stacked X-ray analyses. 

3. We demonstrate using 24 /zm and X-ray stacking 
analyses that galaxies undetected to /24^ m ~ 8 /iJy are 
faint because they have lower SFRs and/or lower ob- 
scuration than 24 /xm detected galaxies, and not be- 
cause they are deficient in PAH emission at a given Ltr- 
We infer that typically half of the bolometric luminos- 
ity of these 24 /xm-undetected galaxies is emitted in the 
UV. Comparing the rest-frame UV composite spectra of 
24 /im-undetected galaxies with those in the top quar- 
tile of detected objects shows that the latter have low- 
ionization interstellar absorption lines that are ~ 2 times 
stronger than the former, indicating some combination of 
more turbulent interstellar media and a larger neutral gas 
covering fraction. This is consistent with our conclusion 
that galaxies undetected at 24 fim have lower SFRs than 
24 /xm-detected galaxies, and therefore a lower input of 
kinetic energy and dust into their ISM. 

4. The 24 /jm (and deep X-ray) data indicate that 
galaxies whose current star-formation episodes are older 
than 100 Myr and have infrared luminosities 10 10 < 
£iR S 10 12 L© appear to follow the local relation be- 
tween rest-frame UV slope and dust obscuration, imply- 
ing that such galaxies at z ~ 2 have moderate amounts 
of dust extinction and that their UV slopes can be used 
to infer their extinction. Galaxies younger than 100 Myr 
have rest-frame UV colors that are redder than expected 
given their inferred £ir, indicating they may obey a 
steeper extinction law. These young galaxies have the 
lowest stellar masses, but span the same range in bolo- 
metric luminosity as galaxies with larger stellar masses. 

5. Galaxies with Lboi ^ 10 12 L Q , including radio- 



detected SMGs, are typically ~ 10 — 100 times more 
dust obscured than their UV spectral slopes would indi- 
cate, assuming their 850 /im-inferred infrared luminosi- 
ties. The 24 /im-inferred infrared luminosities of radio- 
detected SMGs are systematically a factor of 2 — 10 times 
lower than those predicted by their 850 fjxa fluxes; adopt- 
ing the 24 /Ltm estimates implies dust attenuation factors 
that are ~ 5 — 50 times larger than their UV spectral 
slopes would indicate. Regardless, such galaxies will of- 
ten be blue enough to satisfy the U n G7Z criteria, so find- 
ing these dust obscured galaxies in optical surveys is not 
uncommon. 

6. A comparison between the dust obscuration in 
z ~ 2 and z = galaxies suggests that galaxies of a 
given bolometric luminosity are much less dust obscured 
(by a factor of ~ 8 — 10) at high redshift than at the 
present epoch. This result is expected (a) as galaxies 
age and go through successive generations of star forma- 
tion and dust production and (b) if the distribution of 
dust and star formation in galaxies becomes more com- 
pact over time (e.g., through mergers or interactions) 
resulting in greater dust column densities towards star- 
forming regions. We find that star-forming Distant Red 
Galaxies (DRGs) with stellar masses M* > 10 11 M© 
and ages > 2 Gyr have bolometric luminosities and dust- 
obscuration factors similar to those of local ULIRGs, sug- 
gesting that such DRGs, like local ULIRGs, carry rela- 
tively large amounts of dust into their current episodes 
of star formation. 

7. Galaxies with the largest stellar masses at z ~ 2 
(M* > 10 11 M©) also span a large range in bolomet- 
ric luminosity, from those which have red near-IR col- 
ors {{z — A)ab > 3) with little current star forma- 
tion to those ULIRG objects found among optical and 
near-IR selected massive galaxies. Our results suggest 
> 70% of massive galaxies (M* > 10 11 M Q ) in optical 
and near-IR surveys with redshifts 1.5 < z < 2.6 have 
Lboi > 3 x 10 11 L (SFRs > 20 M Q yr" 1 ), comparable 
to and exceeding the luminosity of LIRGs. 

8. Similar to lower redshift studies, we find a trend 
between specific SFR (SFR per unit stellar mass) and 
stellar mass at z ~ 2 which indicates that the observed 
star formation contributes more to the buildup of stellar 
mass in galaxies with low stellar masses than in those 
with larger stellar masses at z ~ 2. This trend be- 
tween specific SFR and stellar mass indicates a strong 
decrease in cold gas fraction as function of stellar mass, 
consistent with results from near-IR spectroscopic ob- 
servations, and suggests that galaxies with large stellar 
masses (M* > 10 11 Mq) at z ~ 2 will quickly cease star 
formation. Combining optical and near-IR selected sam- 
ples, we find a large range in the evolutionary state of 
galaxies at z ~ 2, from those which have just begun to 
form stars and which have large gas fractions to those 
which are old, massive, and have little remaining cold 
gas. 
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